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The field Q(2cos (^)), its Galois group, and length ratios in the 

regular n-gon 

Wolfdieter L a n g ^ 
Abstract 

The normal field extension Q(p(n)), with the algebraic number p{n) := 2 cosf— J, for n e N, 

\ is related to ratios of the lengths between diagonals and the side of a regular n— gon. This has been 

■ considered in a paper by P. Steinbach. These ratios, numbered k = 1, are given by Chebyshev 
polynomials S{k — l,x = p{n)). The product formula for these ratios was found by Steinbach and is 
rc-dcrivcd here from a known formula for the product of Chebyshev S"— polynomials. It is shown that 
it follows also from the S"— polynomial recurrence and certain rules following from the trigonometric 

. nature of the argument x = p{n). The minimal integer polynomial C(n,x) for p(n) is presented, and 

its simple zeros are expressed in the power-basis of Q(/9(n)). Also the positive zeros of the Chebyshev 

■ polynomial S{k — l,p{n)) are rewritten in this basis. The number of positive and negative zeros of 
I C{n,x) is determined. The coefficient C(n,0) is computed for special classes of n values. Theorems 

on C{n, x) in terms of monic integer Chebyshev polynomials of the first kind (called here t) are given. 
i^H I These polynomials can be factorized in terms of the minimal C-polynomials. A conjecture on the 

. discriminant of these polynomials is made. The Galois group is either Zs(^n)j the cyclic group of the 

order given by the degree 6(n) of C{n,x), or it is a direct product of certain cyclic groups. In order 
to determine the cycle structure a novel modular multiplication, called Modd n is introduced. On the 
reduced odd residue system Modd n this furnishes a group which is isomorphic to this Galois group. 
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1 Introduction and Summary 



oo 

O; 

Length ratios between diagonals and the side of a regular n— gon (called diagonal/side ratios, abbreviated 
I DSRs) have been considered by Steinbach [29]. He gave a product formula for these ratios (called by 
him diagonal product formula (DPF)), and for the pentagon and heptagon details were given. In the 
pentagon case the quadratic number field Q(\/5) with the basis < 1,99 > for integers of this field turns 
up. Here if, the golden section, is identified with p(5) which is the length ratio between any of the two 
diagonals and the pentagon side. For the heptagon there are two different diagonal length and the two 
ratios between them and the side length have been called p := p{7) and a. The DPF allowed to reduce 
I all products and quotients of p and a to Q-vectors with the basis < 1, p,a > . For example, = 1 + cr, 
and one can use instead the power basis < 1, p, p^ > of the algebraic number field Q{p{7)) (the use of the 
same acute bracket notation for different bases should not lead to a confusion). The minimal polynomial 
of p(7), in the present work called C(7,x), is — — 2 x + 1, and this was also given in [29]. Therefore, 
[Q(/o(7)) : Q] = 3, which is the degree of this field extension. For odd n, n = 2k + 1, Steinbach gave 
a polynomial with one of its roots p{2k + 1). In the present work the minimal polynomial for p{n), 

called C(n,x), for every n G N, is given in terms of the known one for cos ^"^^ ' called ^{N,x), with 

N = 2n. The connection between ^{N,x) and divisor product representations of A^, and to Chebyshev 
T— polynomials has been worked out earlier by the author. See the links under OEIS [27] A181875 and 
A007955 (in the sequel OEIS Anumbers will appear without repeating this reference). C{n,x) turns out 
to be an integer polynomials of degree 6{n) = A055034(n) . All its zeros are known, and they are simple. 
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Q(p(n)) is the splitting field for the minimal polynomial C{n,x) of p{n). It is a normal extension of 
the rational field. These C-polynomial zeros are written in the power basis of Q(p(n)) with the help of 
certain scaled Cheyshev T— polynomials (called here t with their integer coefficient array A128672 ). This 
array can also be used to rewrite the positive zeros of the Chebyshev S{n — 1, x) polynomial, related to 
the DSRs, in this power basis. The use of 5— polynomial technology allows also for a re-derivation of the 
product formulae (DPF) for the length ratios between the regular n— gon diagonals and side. 

In section 2 S{k — l,p{n)) is shown to yield the DSRs of the regular n-gon. The DPF is also re-derived 
there, and the independent products are extracted. This is called reduced algebra over Q. In section 3 
the minimal polynomial C{n, x) is presented, and its zeros are related to the power basis of the algebraic 
number field Q(/3(n)). In section 4 the Q-automorphisms of this field, the Galois group G{Q{p{n)) /Q) , 
is treated. It is the cyclic group except for infinitely many n- values where it is the direct product 

of cyclic groups, hence non-cyclic. For n from 1..100 there are 30 non-cyclic groups. A novel modular 
multiplication on the odd numbers, called Modd n, is introduced which serves to determine the cycle 
graph structure of this Galois group. The multiplicative group M)dd n is isomorphic to this Galois group. 

2 Regular polygon diagonals/side ratios (DSRs) 

Motivated by a paper of Steinbach [29] we consider the diagonals and the side in the regular n— gon 
(inscribed in the unit circle). The vertices on the unit circle are called V^'^\ A; = 0, ?i — 1, n = 2, 3, 
V^''"^ has coordinates (1,0); the side length (with the radius' length 1 unit) is s{n) = 2 sin (^~^- The 

length of V^^"^ V^f"^ is = 2 sin ^ — ^ , for k = 1, ...,n — 1, hence di = s{n). We only need one side 

and the diagonals in the upper half plane, i.e. it suffices to consider k G {1, [n/2j}. For even n, the 
largest diagonal (k = length 2, lies on the negative real axis. See Fig. 1, and Fig. 2, for the case 

n = 10 (decagon), and for n = 11 (enneagon), respectively. The length ratios of interest, the DSRs, are 

An) 

(n) d, 

R}^'' = -JL^ for the given /e— values. Here Ghehyshev 5"— polynomials in their trigonometric form enter 
df> 

the stage (see A049310 for their coefficients): 

4") = S{k-l,p{n)), (1) 

where we use the second ratio (for the smallest diagonal) p{n) := = 2 cos (^~^- Remember that 

S{n, x) := U{n, x/2) with Gheyshev J7— polynomial (second kind). It is well known that the zeros of the 
polynomials S{n — l,x) are Xi^_^ = ±2 cos ( — j , for k = 1, [^^y^] • Note that in the even n case 

the zero x = appears twice from ±0. The S"— polynomials are orthogonal polynomials defined on the 
real interval [—2,-1-2]. Hence their zeros are guaranteed to be simple. They belong to the Jacohi class. 
Due to (anti)symmetry it is sufficient to consider only the positive zeros (we disregard x = in the even 
n case). Now also Chebysiev T— polynomials (first kind; see the coefficient table A053120) enter, which 
are written in terms of S*— polynomials. (The T— polynomials appear as trace polynomials in a 2 x 2 
transfer matrix when one solves the three term recurrences of these classical orthogonal polynomials.) 

4'^"'^=^&'^ = 2T{k,p{n)/2) =:Up{n)) ^i{k,p{n)) (2) 



Sik,p{n)) - S{k-2,p{n)), k = l,..., 



n 



(3) 



Note that these positive zeros decrease with k, and that tj. is an integer polynomial of degree k because the 
5— polynomials are integer. The coefficient array for these monic i-polynomials is A127672 . In [1] these 
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polynomials are called Chebyshev C-polynomials, but we use the letter C for the minimal polynomials 
of p{n). The zeros can be replaced by DSRs from the above given eq. (1). 



Positive zeros of S(n — l,x) from DSRs: 



X 



(n-l) 



R 



in) 
k+l 



R 



(n) 

k-l ■ 



(4) 





Figure 1: Dekagon n 
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Figure 2: Enneagon (n=ll) 



As an aside we give the factorization of S{n — in terms of the polynomial 
P{k,x) := — 2 cos ( — — ;— with positive zeros. 



2k + l 
S{n-l,x) = x®("-i)p 



n 



1 



,x^ (-i)LVJ p ^ 



n 



1 



(5) 



where 9{n — 1) := if n is odd, and 1 if n is even. This factorization is also considered in [4] for the 
U{n, x) polynomials for even and odd 7i separately. For example, for the heptagon n = 7 one finds the 
factorization of S{6,x) with the following real polynomial P{3,x). 



P{3,x) 



(2a(7) -l)x^ + 2p{7)x - 1 



(6) 



with p{7) = i?2 '^C'') ■— ^3' also the introduction). The three positive zeros are, written in 

terms of the DSRs, p{7), a{7) — 1 and cr(7) — p{7). Here we used the general translation formulae eq. 
(4) to rewrite the positive zeros of S{n — l,x) in terms of DSRs. One should also write this in terms of 
powers of p{7) with the help of eq. (1) (see eq. (7) below). 
n — 1 \ 

— - — ,xj can be considered as characteristic polynomial for a ( [^^y^J + l)-term recurrence 

sequence {fk}'kLo integral domain of Q(p(n)) (represented by a 5(n)-tuple of ordinary (rational) 

integers (Ai^k, -^2,k, As{n),k)- Here S{n) is the degree of the minimal polynomial of p{n) which will be 
discussed in sect. 3. In fact, one has |_^^J such (5(n)-tuple sequences corresponding to the independent 
inputs. As a simple example take the n = 5 pair ((5(5) = 2) of sequences corresponding to P{2, x) = 
— (2 (/9 — 1) x + 1, with the golden section 99 := p{5) and the simplest input. This is S{k, ^/b) = Ai^k 1 + 
^2,fc with input 5(-l, ^/S) = 0, 5(0, y/h) = 1. See A005013 (fc + 1) (-1)'=, and 2 A147600 (fe - 1), for 
Ai^k and ^2,^, respectively. 

For later purposes it is also useful to give a dictionary between the positive zeros of S{n — 1, x), i.e. those 
of -P(L^^J ,ic)j and powers of p{n). One direction follows from the above given formula, eq. (2), for 



Each P 
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(n-1) 



in terms of the monic integer polynomials i{k,p) = 2T{k, p/2). The dependence on n has here 



X 



been omitted. 



Positive zeros of S(n— l,x) from p(n)-powers: 
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p{n) 



k-2l 



(7) 
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(omitting n). For 



The coefficients are given in the triangle A127672 . E.g., = 5p — 5 p^ 
the heptagon this shows that the three positive zeros of 5(6, x) are p(7), (t(7) - 1 = p2(7) _ 2 and 
(t(7) — p{7) = /o(7)^ — p(7) — 1 . In this case these three zeros could be used as a Q-vector space 
basis. Note that the n dependence is via the i-polynomial variable p{n) (which satisfies C{n, p{n)) = by 
definition of the minimal polynomial C). For example, it is true for all n that the second largest positive 



zero of S{n 



- 1, x), namely 



(n-l) 



, is always p{n) 



2. Hence, from above, R 



(n) 



1 = /o(n)^ - 2; or 

if one calls R'^'' =: cr{n) then p{n)^ = (j{n) + 1. This has been noted for n = 7 above but is holds 
in general, showing that o"(n) is algebraically dependent on p{n). This fact has been observed already in 
[29] where it appeared as a special product formula (see the later DPF eq. (13), m = k = 2, or eq. (11) 
with k = 2). 

The inverse formulae (DSRs in terms of p-powers) have already been given above in eq. (1). Here one 



uses the S'-triangle A049310 for the translation. E.g., R!^ 



(n) 



-2p{n) + p^{n) = p{n){p{n) 



2) 



= p{n) {(y{n) — 1) . Of course, R)l''' is interesting only for n > 9 because only then the corresponding 
diagonal lies in the upper half-plane (including the negative real axis). 

Note that the number of positive zeros of S{n — \,x), i.e., |_^^J , is less or equal to |_^J, the number 
of DSRs for diagonals in the upper half plane including the negative real axis. The zero x = 1 never 
appears. As can be seen in the n = 7 case these zeros can nevertheless be used as Q-vector space basis (in 
sect. 3 it will become clear that the degree (5(7) is also 3). In the case n = 9 the two numbers are both 4, 
but the degree (5(9) is 3, hence only three of the zeros and three of the DSRs are rationally independent. 



To wit: x(^) 



5xf^ + 3x\"' - Ax'^' and R, 



.(8) 



(8) 



,(9) 



R 



(9) 



+ rP 



1 + p{9). 



Because of symmetry only DSRs for diagonals of the upper half plane (including the negative real axis 
n 



t.e. 



k G {0,1,..., 



accomplished by the ru 



} , are of interest. The reduction for other k values, also negative ones, is 
es 



-R 



in) 



^-\k\ 



-R 



in) 
k\ 



ii) for A; G { - +1, n} 



R 



in) 



(8) 



These rules follow from eq. (1) with the trigonometric definition (with the specific value of p{n)) of the 
^-polynomials, and therefore also negative values for the DSRs show up (the interesting DSRs are, of 
course, positive), i) follows also from S{—\n\,x) = —S{\n\ — 2,x) which derives from a backward use of 
the recurrence (given later in eq. (12)). This rule does therefore not depend on the special choice of the 
variable x, in contrast to the rules o) and ii). With negative k one counts the diagonals in the clockwise 
direction, ii) is used to translate from the lower to the upper half-plane. 

The product formula for Chebyshev S'-polynomials of different degree, but with the same argument, is 
well known. See e.g., [1], p. 782, 22.7.25, for the [/-polynomials, and replace the T-polynomials by the 
S'-polynomials via the trace formula given in eqs (2) and (3). 



S{m — 1, x) S{n — 1, x) 



1 



4((f)2 



[S{n + 'm,x) — S{n + m — 2,x) — 

S{n — m,x) + S{n — m — 2,x)] n > m . (9) 



For our purpose x = p{n) = R^'^ ^ ±2 but its special form is for this formula not of interest . It is in 
fact symmetric under exchange of n with m if the rules for negative indices on the S'-polynomials, stated 
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above, are employed. Here we restrict to n > m in order to have non-negative indices. Thus the product 
formula for the DSRs, now with x = p{n), but without using its specific value, is 

(4 - = - 4-^^-1 - <Ui + ,k>m>0. (10) 

For m = this becomes trivial. This is not yet the DPF given by Steinbach in [29] which linearizes 
the product of Rm^ ^^k^ ■ order to eliminate the pre-factor (4 — p{n)'^) one can use the three term 
recurrence relation of the orthogonal 5-polynomials, written for the DSRs by eq. (1) as a special product 
formula 

[Rec,k](") : = ijg, + , ^ > . (11) 

Prom now on we also use the specific form of p{n). Thus also rules o) and ii) of eq. (8) will be 
applicable. For k = one uses (see ii)) R^_^l = — ii^"^ = —1, and it becomes trivial. For k = 2 
it shows that p{n)'^ = cr(n) + 1 for all n, if one uses cr(n) := ■ This recurrence can now be 
used twice as R^2^^ (-^2"^ -^i"^) ™ ^l- (10) to produce the following recurrence for two step differences of 

{pk+2 - Pk) - {Pk - Pk-2) = Ck - Ck-2 , (12) 

where we used the abbreviation Ck = c^"'"*^ := R^^J^_^-^ — R^i^].^_^i- This shows that pk+2 — Pk — Ck 
is k independent. From the inputs P2 — Po — co = 0, due to the recurrence relation eq. (11), and 
Pi — P-i — c^i = with the help of the rule i) from eq. (8), this leads to the recurrence Pk+2 — Pk = Cfc, 

(n) 

with the inputs po = and pi = —pi = Rm ■ The solution for pk is found for even and odd k separately, 
where again the rules for negative indices are employed. Both solutions can then be combined as 

DPF [m,k]H : = <V(2,+i) ,l<m<k. (13) 

j=0 

This is finally the DPF found by Steinbach in [29] and gives a linearization of the DSR products. We 
need only to consider m > 2 {m = 1 becomes trivial), and due to the symmetry of this formula under 
the transformation k ^ n — k (using the rules of eq. (8)) is suffices to consider k G {1, 2, ... }. Hence 

f\-\\ 

one has only to consider the ( P^'O'^^'^^^ 2 < m < k < [^J . 

The DPF formula looks un-symmetric with respect to m -f-)- k, but it is, in fact, symmetric because 

m— 1 

R^j^^^_(^2j+i) for A; < m reduces, due to cancellations after using the rule i) from eq (8), to the 

j=0 

expected sum with only k terms. One can see this for even and odd m > k separately, remembering 
that iJg""^ = in the former case. 

The idea is to work out, for a given n > 4, all DPFs of interest, using the rules from eq. (8), especially 

i) and ii) in order to write all ( 2 ) P'^oducts as linear combinations of the DSRs. 

Example 1: n = 7 heptagon, treated in detail in reference [29]. Here we recapitulate and link to OEIS 
[27] sequences (the analoga of Fibonacci numbers in the pentagon case). Superscripts and arguments 7 
are suppressed. [2,2]: p"^ = a + 1. [2,3]: pa = a + p + 0. [3,3]: a"^ = p + a + 1. Here, R5 = R2 
was used. This shows that the Q- vector space basis is at most < 1, p,a >. It will be shown to be indeed 
the heptagon basis in sect. 4, where the power basis will be used instead. 

With these DPFs one can compute all powers of interest in the heptagon basis < l,/5, o" > (this has 
already been done explicitely for a in [29], p. 28). 

= A052547 (fc-2)1+ A052547 (fc - 1) p+ A006053 (fc) a. A: > 0, 
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p-k = A077998 (fc) 1+ A077998 ffc- 11/9- A006054 (A: + IW, > , 

= A106803 fA: - 1) 1+ A006054 (A; - 1) /o + A106803 (/c) cr. fc > 0, 
o--^ = {a - pf = A052547 (fc) 1 - A006053 (/c + l)p- A052547 (fc - 1) cr, A: > 0, 
{paf = {p + erf = A120757 (fc) 1+ — A006054 (fc - 1)1 p + 4 Al81879(fc) a, /c > 0, 
(go-)-fc = + = A085810 (fc)l + (-1)^ A181880 (fc-2)p + f-D^+i A116423 (fc + 1) a, fc > 0. 

One can also compute = ^1+ B p + a and find with N{a, b, c) = — 6^ — — 

al + 6p + cp 

2a 6^- a6c + 0^6 + 6^0 + 2a'^c - ac^ + 2b 

^ = M( \ ^ - 62 + + 2ac - 6c) , B = \ . (6^ - + a6) , (14) 
N{a,b,c) N{a,b,c) 

^ = AT( \ \ - 6^ + ac - 6c). (15) 
A* (a, 6, cj 

Example 2: n = 9 (Enneagon), where we use p = a = R^3\t = R^^'^- [2,2] : p"^ = 1 + a, 

[2,3] : pa = p + T, [2, 4] : pr = a + t, [3,3] : cr^ = 1 + r + ct, [3, 4] : ar = p + a + r, 
[4, 4] : = 1 + p + 0" + r. This DSR-algebra (over Q) shows that not all ratios (including = 1) 
are linear independent: t (p — 1) = a = p^ — 1 = {p + 1) (p — 1), i.e., t = p + 1. Therefore, 
p'^ = p{a + 1) = 2 p + T + 2 = 3 p + 1. In sect. 3 we will see that the algebraic number p{Q) 
has degree 3, and its minimal polynomial is indeed C(9, x) = — 3x — 1. The enneagon basis is 
thus < l,p, o" > which can be related to the power basis < l,p,p'^ > (remember that we use the same 
notation for different bases). The reduced DSR-algebra (the algebra modulo C(9, p) = ) is = 1 -|- o", 
po" = 1 + 2p, and a'^ = 2 + p + a. 

Remark 1: [2,2] in eq. (13) becomes {R'^^f = R^'^ = 1, showing that for all interesting values n > 4 
one has (7{n) := R^^ = p{n)'^ — 1 . This is known already from eq (1) for k = 3, and it will be used in 
sect. 3 as iJg"^ rewritten in terms of the power basis of Q{p{n)). 

Remark 2: For n = 4 and n = 5 the second diagonal (defining R^^^) is not of interest because it is 

in the lower half plane. According to eq. (8), ii) i?2 — -^3 ) hence cr(5) = p(5) and the general 
relation [2, 2] from eq. (13) between (j{n) and p(n) leads to p(5) = p(5)2 — 1, the golden section formula 
p(5) = (p := (1 + ^/5)/2. In sect. 2 it will be seen that the minimal polynomial for p(5) is indeed 
C(5,x) = x"^ - X - 1. 

Instead of deriving the DPF from the product formula for S-polynomials and the rules of eq. (8) one 
can prove it directly by induction from the S"— recurrence including negative indices in accordance with 
these rules. 

Proposition 1: 

The DPFs [m,/c]("), eq. (13) but now for all 2 < m, 2 < k, follow from the recurrence, eq. (11), and 
the rules i) and ii) from eq. (8). 

Proof: This is shown by double induction. First one shows this for given m > 2 for all k > 2 by 
induction over k. Then by induction over m for all A; > 2. 

The /c-induction uses as starter [2,2]^"^ which is the recurrence. Assume that [2,A;']*-"'^ is true for all 
k' = 2, 3, k — 1. Now (we omit the superscripts) R2 Rk = R2 {R2 Rk-i — Rk-2) from the recurrence. 
Then use for both terms the induction hypothesis. Note that one obtains for the first term k — 1, and 
for the second one k — 2 terms. However, due to rule ii) only two terms survive in each case. (The same 
type of cancellation was at work when we remarked above on the symmetry of the DPF formula in m 
and k.) Therefore, one obtains R2 {Rk — Rk-2 + 0) — {Rk-i — Rk-3 + 0) which becomes, after use of 
the recurrence applied twice, Rk+i + Rk-i- This is indeed the desired result for [2,k] if one uses again 
rule ii) to get a truncation of the sum after two terms. 

The m-induction uses as starter [2, k] for all k > 2, which has just been established. Then assume that 
[m',k] is true for all m' = 2,3, ...,m — 1. RmRk = R2 Rm-i Rk — Rm-2 Rk from the recurrence. 
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Assume the induction hypothesis for each term, obtaining, after use of the recurrence for the first sum, 

fc-l k-l 

(i?m+A:-(2j+i) + -^m+fe-2-(2j+i) ) ~ (-^m-2+fc-(2j+i) which is indeed the assertion after cancellation 

j=0 j=0 

of the last two sums. □ 
As mentioned in the example 2 the DSR-algebra turns sometimes out to be reducible because some DSRs 
can be expressed as rational linear combinations of other ones. Later it will become clear that this 
happens precisely whenever [^J — 6{n) > 0, and this is the number of linear dependent DSRs. See 
Table 1 for details for n = 3, 12. 

We shall see in sects. 3 and 4 that it is simpler to use the power basis of Q_{p{n)) and the minimal C{n, x)- 
polynomials of the algebraic number p{n) instead of the DSR-algebra. Then one obtains automatically 
the reduced algebra. 



3 Minimal polynomial of p(n), its zeros, absolute term and factoriza- 
tion 

The minimal polynomial of an algebraic number a of degree da is the monic, minimal degree 
rational polynomial which has as root, or as one of its roots, a. This degree da is 1 iff a is rational, and 
the minimal polynomial in this case is p{x) = x — a. For the notion 'minimal polynomial of an algebraic 
number' see, e.g., [22], p. 28 or [25] p. 13. 

For the algebraic number p{n) := 2 cos (—\ for n G N, the degree is 6{1) = 1, and 6{n) — ^ 



_ ^ -7 ~\ / o 

for n > 2 , with Euier's totient function ip{n) = ylOOOOlO (n). This is the sequence A055034 (n). The 
sequence of allowed 6 values is given in A207333 . The array with the indices of the polynomials for given 
allowed 5 values is shown in A207334 . Of course, 6 is not multiplicative, e.g., 2 = J(6) ^ 5(3) (5(3) = 

1-1 = 1. These polynomials can be obtained from the ones of cos ( — ) which are found, e.g., under 

\n J 

A181875 / A181876 . and they have been called there ^'(n,x). See also [17], and [22], Theorem 3.9, p. 37, 

for the degree d(n) of these polynomials. From the trivial identity cos ( — ) = cos ( — ) one finds the 

\n/ \2n 

minimal polynomial of 2 cos called here C{n,x), from 



Cin,x) = 2''(")^' (2n,|j . (16) 

Therefore the above formula for 5{n) derives from d{2n). These polynomials are given for n = 1, 30 
in Table 2, and the first 15 rows of their coefficient array are shown in Table 3. Concerning the parity of 
the degree 6 one has the following lemma. 

Lemma 1: Parity of the degree S 

6(n) is odd iff n = 1, 2, and p^, with an odd prime p of the form 4k + 3, k G Nq, and e > 1. 

Proof: The case n = 1 is clear by definition. For n = 2 one uses the prime number factorization n = 

N 

'^^ nj^i 'P^j -I with /3 > and = 1. It follows from the definition that b{n) = 2^~^ [pj — l)Pj^ 

i=i 

For P >2 this is always even. For /3 = one needs = 1, otherwise this will be even. For A^ = 1 one 
p — \ 

needs — ^ — to be odd, i.e., p = 'ik + 3 with k E Nq- If /3 = 1 one needs A^ = in order to have an 
odd value. Thus n = 2. □ 
In order to relate to sect. 2 we are interested here in n > 4. It turns out that these polynomials are in 
fact integer (not only rational) polynomials. The proof will use the following lemma based on our divisor 
product representation paper [14]. 
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Lemma 2: Minimal polynomial C(n,x) written as a rational function. 



q{dq{n),x) 

with monic integer polynomials p and q with a certain degree dp{n) and dq{n) = 5p{n) — 6{n), respec- 
tively. 

Proof: ^ (2n, I) is obtained from the unique divisor product representation dpr{2n) defined in [14] 
by replacing each a{k) in the numerator, as well as in the denominator, by t{k, |) which is given as a 
difference of monic integer polynomials t which have been defined already in eq. (2), multiplied with a 
certain prefactor. 



27 



+1 



k- 
2^ 



+ 1,2;) - t(| - 1,2;)) if A; is even , 

(18) 

- (t(^,x) - i{^,x)) if A; is odd. 



The monic ^(2n, x) polynomials have degree d{2n) (see e.g., [13]) which implies that all the prefactors 
in this replacement of the a(fc)s in dpr{2n) have to become 1/2'^^^"). (This could be formulated as a 
separate lemma). E.g., n = 34 with dpr{34:) = ^1^7) 0(2) from the numerator the factor 1/2^^"'"^ 
and from the denominator 2^+^, fitting with 1/2'^^'^^) because (i(34) = 8. Therefore, one may in the 
calculation of C(n, x), found above from ^(2n, |), forget about these prefactors in the replacements 
altogether. Thus the numerator, resp. denominator, is a product of monic integer polynomials t leading 
to the monic integer polynomials p, resp. q, of a certain degree dp{n), resp. dq{n). (One could give 
more details on these degrees but this is not important here. Trivially, dp{n) — dq{n) = d{2n). For 
the given example n = 34: (18 + 1) — (9 + 2) = 8 from the degrees of the t polynomials.) Because 
C(n,x) is a minimal polynomial this rational function allows polynomial division without remainder. 
□ 

Proposition 2: C G Z[x] 

C{n,x), the minimal polynomial of p{n) = 2 cos is an integer monic polynomial. 

Proof: From lemma 2 we have q{dq{n),x) C{n,x) = p{dp(n),x) which leads by induction to the re- 
sult that the integer coefficients of the monic polynomials p and q imply integer coefficients for the 
monic polynomial C. Call these monic polynomials q{N,x) := X^^q 9/x', C(M, x) := X]fcl=o '^k^^ and 
p{N + M,x) := X^^o*^ Pi^'' with qat = 1, cj\/ = 1 andpAr+A/ = 1. Collecting terms for x^"*"^^"-^ , for j = 
(0), 1, 2, M, one obtains the formula for the C-coefficient cm-j in terms of lower indexed ones: cm-j = 
j 

PN+M-j — CM-j-k QN~k- With this the inductive proof on j, using the integer coefficients of q,p and 

k=l 

the integer higher C-coefficients due to the induction hypothesis becomes obvious. The starting point is 
the trivial j = case. 

Next we give all the zeros of C(n, x). This follows directly from the knowledge of all zeros of ^{2n, |) 
(see, e.g., [13], from p. 473 of [30]) 



Proposition 3: All zeros of C 

C{n,x) = Yl (x-2cos(^)),nGN. (19) 



" ^ ' ^ TT k 



gcd(k,2 n) = l 



In the product the index k starts with 1 because gcd{0,2n) = 2n, and it stops at n — 1 for all n > 2 
because gcd{n,2n) = n. One has to omit the even k values > which leads to (for n = 1 one takes 
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only / = in the following product) 



C{n,x) = II ("x - 2 cos ("^^i^itll") ) . (20) 



1^0 ^ ^ 

gcd(2 l + l,n) = l 

Proof: This is a simple consequence of the mentioned known results for ^(2 n, |). It implies, en passant, 
a formula for the known degree 6{n) (see A055034 ) in terms of the number of factors of this product. 
□ 

The following proposition lists the nonnegative and negative zeros of C for prime n. A vanishing zero 
(which will be counted later as positive) appears only for n = 2. 

Proposition 4: Non-negative and negative zeros of C(p,x) 

i) If n = 2 then C(2, 0) = , and this is the only case with a vanishing zero. 

ii) If n is an odd prime l(mo(i4) then the number of positive and negative zeros coincides, and this 
number is These zeros are (-1)^+^ 2 cos (^), with k values 1, 2, = 6{n). 

iii) If n is an odd prime 3 (mod 4) then the number of positive zeros exceeds the number of negative 
ones by one, and the number of negative ones is These zeros are (—1)'^'^^ 2 cos (^) with k values 
1,2, and an extra positive zero appears for k = = 6{n) . 

Proof: 

i) Vanishing zeros require n = 2 (mod 4) from ^ = \ and k odd in eq. (20). A vanishing zero appears 
for n = 2. For other such n values, n = AK + 2, K > 1, the odd k = ^ > 1 divides n, hence it does 
not appear in the product of eq. (20). 

p-1 

For odd primes n = p = 2q + 1 the product in eq. (20) is unrestricted, and there are (q—l) + 1 = — - — 

factors, in accordance with the degree 6{p) = The cos function will produce negative zeros for 

2/ + 1 1 

> -, i.e., for 2/ > q. 

p 2 

n)n = p = 4:K+l means that I = 0, 1, K — 1 lead to positive zeros, and / = K, K-\-l, 2K — 1 

Tl — 1 

to negative ones. In both cases there are K = — - — such zeros. For the / values leading to negative 

zeros one can use the formula cos ( — (2 / + 1) ) = — cos i — (p — (2 / + 1)) ) . In this way even multiples 

\p J \P J 

TT 

of — appear, and one produces, counted backwards, beginning with the largest / value, the new values 
P 

2,A, ...,2 K . This proves ii). As a test consider n = 13, X = 3: the 21 + 1 values for the positive 
zeros are 1, 3, 5, and the ones for the negative ones are 7, 9, 11. The latter become (we use underlining 
to indicate that they come with a minus sign in the final formula) 6, 4, 2. Rearranged these values are 
1) 2, 3, 4, 5, 6, as given in ii). 

iii) For n = p = AK + 3 a similar analysis leads to ^ = 0, 1, K, and I = K + \, 2K for positive, 
and negative zeros, respectively. This shows that the number of negative zeros is = — ^ — , and the 
number of positive ones exceeds the one for negative zeros by one. Again, the odd values leading to nega- 
tive zeros are transformed to even ones 2, 4, 2K with a minus sign in front of 2 cos. Take, e.g., n = 19, 
K = A, with the new k values of iii) 1,2, 3, 4, 5,6,7,8, 9. □ 

For general values n one can also find the number of positive and negative zeros of C{n,x). In the 
non-prime or non-power-of-2 case the gcd restriction in the product excludes certain / values. Therefore 
one has to eliminate from the unrestricted 21 + 1 values in the product of eq. (20) all odd multiples 
for each odd prime dividing n. Multiples of 2 are irrelevant for even n , therefore only the odd primes 
dividing n are of interest. This shows that the prime factors of the squarefree kernel of n, denoted here by 
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sqfk{n), i.e., the largest squarefree number dividing n (see A007947) will be of interest. This squarefree 
kernel is also known as radical of n, denoted by rad{n). We will denote the set of primes of this kernel 
by sqfkset{n), E.g., sqfkset{2^ • 3^ • 11) = {2,3, 11}, sqfk{2^ ■ 3^ • 11) = 66 (strip off all exponents in 
the prime number factorization of n). Because one may encounter multiple counting {e.g., 15 is hit by 
the multiples of 3 as well as 5 for any n which has in its squarefree kernel set the primes 3 and 5, and 
which is larger than 16) one can employ PIE, the principle of inclusion-exclusion, e.g., [5], p. 134, to get 
a correct counting. 

Proposition 5: Number of positive and negative zeros of C(n, x) 

i) If n = 2"*, m > 1, one has for m = 1 a vanishing zero (see proposition 4 '>-))■ For m > 2 the number 

Tl 

of positive zeros, which is — , coincides with the one for negative ones. This is in accordance with the 



n 



/vr 



degree 5(2"") = 2""-^ = ^. These zeros lie symmetric: ± 2 cos ( - (2 Z + 1) ) , / = 0, 1, ... 



n 



1. 



2 " \n 

ii) If n = 1 then there is only the negative zero —2. If n > 3 is not a power of 2 then the number of 
positive zeros, called S+{n), is 



M{n) 

in) = K{n) + i-^r 

T = l 



E 

<il,i2,...,ir> 



Lin) 



Ph--- Pi 



(21) 



where M{n) (sometimes called uj{n)) is in the odd n case the number of elements (cardinality) of the 
set sqfkset{n), denoted by \sqfkset{n)\, (see A001221 ). M{n) is for even n the number of odd primes 
in sqfkset{n), i.e., \sqfkset{n)\ — 1. This is because multiples of 2 are irrelevant here. The sum 

< ii,i2, ...,ir > extends over the ( ) combinations 1 < zi < ^2 < ... < ir ^ M. Only the odd 



primes from the set sqfkset{n) enter. The values K{n) and L{n) depend on the parity of n and they are 
given by 



a) neven : K{n) ■ 

131) n odd, \{mod\) : K{n) 
/32) n odd, 3{mod4) : K{n) 



n 



4 

n — 5 
4 

n — 3 



L{n) 

L{n) 
L{n) 



2K{n) + 1 

2K{n) + 1 
2K{n) + 1 



n 



2 

n — 1 



(22) 

(23) 
(24) 



Note that also negative values for the floor function may appear. In this way the pure prime case from 
proposition 4 is also included. 

iii) If n > 3 is not a power of 2 then the number of negative zeros, called 6-{n), is 



M(n) 

N{n) + ^ (-1)^ 

r=l 



5.{n) 

where M[n) and L{n) are as above, and 



<il,i2,...,ir> 



P{n) 

Ph ■■■Pi. 



L{n) 



Pii ■■■Pi 



a] n even : 



N{n) 

pi) n odd, l{mod4) : N{n) 
132) n odd, 2,{inodA) : N{n) 



n 



2 

n — 1 
4 

K(n) 



n 



n 



4 

P{n) 
- 3 



Pin) 



n 



n 



P{n) 



n 



2. 



(25) 



(26) 

(27) 
(28) 



Proof: 
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n 



-\m—l 



1 values contribute to the product of eq. (20), compatible with the 



n 



. The negative zeros appear for 2 Z + 1 = 2 2™ + 1 



, n- 



1, hence 



i) Ah I = 0, 

degree <5(2™) ^ _ ^ 

there are ^ of them. With the cos formula given in the proof of proposition 4 ii) with p — )• n, they can 
be rewritten in the notation, also used in the above context, when read backwards, as 1^,3, 2™~^ — 1 . 
Therefore one obtains the same cos arguments like for the j positive zeros, only the overall sign is 
different. 

ii) For n > 3, not a power of 2, one uses PIE to count the positive zeros. For this consider the odd 
multiples of some odd number a, say (2fc + 1) a, (only such multiples come up as 2/ + 1 values in 
the product (20)) up to k^ax such that all positive zeros are reached. If the largest 21 + 1 value in the 



product which leads to a positive zero is n, then k„ 



n 



and there are kmax + 1 of these 



odd multiples of a. In the following application of PIE a will be taken as any odd prime or product of 
odd primes from the set sqfkset{n). 

a) Case even n, not a power of 2: The counting of the positive zeros starts at level r = with all 

n- 2 



positive ones in the unrestricted product eq. (20). There are 
n - 2 



+ 1 of them (the I values are 



0, 1, 



exceeding 2 



n 



). In the next step, r 
- 2 



This leads to a subtraction of the form 



1, all odd multiples of every odd prime pi^ in sqfkset{n) not 
+ 1 are discarded (there are M{n) such primes, where M{n) := \sqfkset{n)\ — 1). 

L{n) 



M{n) 

E 

n=i 



1 



+ 1 with L{n) given in the proposition 



ii) a). Now double subtractions may have appeared and one adds, at step r = 2, all odd multiples of 
the product of two odd primes from sqfkset{n), call them Pi^Pi2 with ii < 12 (because no problem of 



«2) 



This leads to the term + 

l<ii<i2 <M{n) 

Now one has to correct for triple products in the step r = 3 , etc., up to level r 



L(n) 



1 



+ 1. 



over-subtraction appeared for ii 

Pii Pi2 

-- M{n). This leads 

to the formula given in part ii) a) if one observes that the +1 term in each r-sum becomes a {^^^^ 
term because of the number of terms of this sum. Then, because the alternating sum over row M{n) in 
PascaFs triangle (see A007318 ) vanishes ((1 — l)-'^(") = 0), all these terms add up to —1 (the missing 

n — 2 

r = term). Therefore the formula starts with — - — , the K(n) value given in ii) a), and only the 
floor term remains for each r-sum . 

(51) The counting of the positive zeros in the odd n case distinguishes the two cases l(mo(i4) and 
3 {mod 4). Here we discuss the former one. The procedure is the same as the one given for case a). 
One has just to determine the boundary values for the Imax value leading still to a positive zero. This is 

Imax = K — 1, if n = 4K + 1, i.e., — - — . Hence there are Imax + 1 such / values to start with at level 

r = 0. Again the extra +1 can later be taken as r = term ("^'^q^^) for the vanishing alternating sum 
over row No. M(n) in Pascal's triangle. Here M(n) := \sqfkset{n)\. Therefore the PIE formula starts 



n 



+ 1 



n 



with the K(n) given in f31) of the proposition. The L{n) of the PIE sum is now 2 ^ , _ ^ 
the largest 2 / -|- 1 value leading to a positive zero. 

(32) In the case M{n) := \sqfkset(n)\, n = AK + 3 the maximal 21 + 1 value leading to a positive 

71 — 1 

zero, the L{n) in the formula, \s2K + 1 = , coming from the largest / value which is K. This K 

is the K{n) in the formula claimed for this case /32). 

iii) For the number of negative zeros the counting is done by finding all the odd multiples of odd primes, 
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or products of them, which cover the 21+1 range for the negative zeros. This is the difference of the 
number of such multiples for the whole range and the range for the positive zeros. 

a) If n is even, not a power of 2, the whole range is determined by Imax + 1 = n — 1. From the above 
ii) a) case one knows the corresponding upper bound for the positive zeros, therefore the PIE formula 



77 — 2 

has N{n) = (^1-^ + 1) - ( 



n 



+ 1), the number of factors in the unrestricted product which 



lead to negative zeros, which is the value claimed in the proposition iii) a). Accordingly, P{n) = n — 1 
and L{n) = 2 |_^J + 1 as given in eq. (26). 

In this n = l{mod4:) case the maximal number 21 + 1 in the product is n — 2, determining 
P{n). The corresponding number L{n) for the positive zeros is taken from above as — ^ — . There are 

N{n) = ^ — ^ — + 1^ — ^ — ^ — + 1^ = — ^ — unrestricted factors in the product with negative 
zeros. 

j32) In this n = 2>{mod4) case one has also P{n) = n — 2 like for any odd n. L{n) is taken from the 

/ 77, 3 \ I Ti 3 \ 77 3 

corresponding ii) case as — ^ — , and N{n) is obtained from ( — ^ — + 1 j — ( — ^ — + 1 j = — ^ — . 
□ 

We give the first entries of the 6+{n) and (5-(n) sequences. Remember that the vanishing zero for n = 2 
is here counted as positive. The A-numbers given for 5-{2k) and 5-{2k) are conjectured. 
n even: 

6+{2k), k = 1,2,...; A055034 : [1,1,1,2,2,2,3,4,3,4,5,4,6,6,4,8,8,6,9,8,6,10,...], 

5-{2k), k = 1,2,...; A055034 with first entry : [0,1,1,2,2,2,3,4,3,4,5,4,6,6,4,8,8,6,9,8,6,10,...], 

n odd, l{mod4:): 

6+{4k + 1), A; = 0,1,...; [1,5,9,13,17,21,25,29,33,37,41,45,49,53,57,61,65,69,73,...], 
(5_(4A; + 1), A; = 0,1,...; [1,1,2,3,4,4,5,7,6,9,10,6,11,13,10,15,12,12,18,16,14,...], 
n odd, 3{mod4:): 

6+{ik + 3), A; = 0,1,...; [1,2,3,2,5,6,5,8,6,6,11,12,8,10,15,10,17,18,10,20,...], 
5-{4k + 3), A; = 0,1,...; [0,1,2,2,4,5,4,7,6,6,10,11,8,10,14,8,16,17,10,19,20,...]. 
The start of the sequences for odd n is therefore 

5+{2k + 1), k = 0,1,...; [0,1,1,2,1,3,3,2,4,5,2,6,5,5,7,8,4,6,9,6,10,11,...], 

5- {2k + 1), A; = 0,1,...; [1,0,1,1,2,2,3,2,4,4,4,5,5,4,7,7,6,6,9,6,10,10,6,11,...]. 
Finally, combining for all n: 

6+{n), n = 1,2,...; A193376 : [0,1,1,1,1,1,2,2,1,2,3,2,3,3,2,4,4,3,5,4,2,5,6,4,5,...], 

6- {n), n = 1,2,...; A193377 : [1,0,0,1,1,1,1,2,2,2,2,2,3,3,2,4,4,3,4,4,4,5,5,4,5,6,...]. 
On can check that 6{n) — {6+{n) + 6-{n)) vanishes in each case. 

Observe that from these examples it seems that 6^{2k) = 5-{2k) =A055034(A;) = 6{k) for k > 2. 

ipUk) 

This is compatible with 5{2k) = 26{k), i.e., — - — = 2ip{2k) for these k values. The latter eq., which 
holds also for k = 1, can be checked by considering the two cases k = 2'^^^^ ■ (odd number) and k odd. 
Recall that 6{2) = 6{1) = 1, because d{l) was special (it is not ^ ). 

Zeros of C(n, x) in the povi^er basis < p(n)° = 1, p(n)^, p(n)^(")~^ > 

As will be explained in the next section, the minimal polynomial C{n,x) with degree 5{n) leads to the 
splitting field Q{p{n)) which is a simple field extension of the rational field Q, and the degree of Q{p{n)) 
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over Q, the dimension of Q{p{n)), considered as a vector space over Q, is just the degree of C. In standard 
notation [Q(/3(n)) : Q] = 6{n). Therefore it is interesting to write the zeros of C in the power basis 

for this Q-vector space. This task is accomplished by using for the zeros |W = cos (-(2/ + 1)), / e 



n 



{0,1,..., 



n 



2 



} with gcd{2l + l,n) = 1 (see eq. (20) ), the formula ^{"^ = i{2l + l,p{n)) = 2T{21 + 



p{ti\ 

1, ) (compare with eq. (2)). Remember that the integer coefficient array for these t-polynomials is 

shown in A127672 . Of course, one has to reduce this integer polynomials using C(n, p(n)) = 0, i.e., one 
has to replace p{n)^^^^ (and higher powers) by an integer polynomial of lower degree. This can be done on 
a computer and Table 4 has been found with the help of MaplelS [18]. For example, the zeros of C(15, x) 
are^T^ = p(15), |f ^ = -2 + 3 p{15) + p(15)2 - p(15)3, |f ^ = -1 - 4p(15) + p(15)3, |f ^ = 
2 — /3(15)^ . Here the reduction has been performed with /o(15)^ — t- — p(15)^ + 4p(15)^ + 4p(15) — 1. 
One can check these zeros which have approximate values 1.956295201, 0.209056928, —1.338261216 and 
— 1.827090913, respectively. We will call the zeros with increasing value ^1,^2, •••i^5(n)- 

Our next objective is to compute C{n, 0) (the sign and the number multiplying x", also called the absolute 
term) and relate it first to the cyclotomic polynomial cy{n, x) (the minimal polynomial of the complex 
algebraic number , see e.g., [30] or [9], p. 149, Exercise 50 a and b), evaluated at a; = —1, which is 

n.-l 

cy(n,-l) = (-1)'^(") n (1 + ^'"'") > ^ = 2, 3, ... (29) 

gcd{k,n) — l 

We used the fact that the product in the definition of cy{n,x) has ^{n) factors (the degree as minimal 
polynomial). For n = 1 one has cy{l,—l) = —2, which fits this formula if one defines := 1 and 

takes the undefined product as 1. One may rewrite this for n > 2 by extracting e~ , using for the sum 
of the (p{n) terms the formula 

4 E ^ = ^ , n = 2, 3, ... (30) 

fc=i 

gcd(k ,n) — l 

This is a formula listed by R. Zumkeller under A023022 which is found in[2], p. 48, exercises 15 and 16, 
written such that both sides depend only on the distinct primes in the prime number factorization of n. 

99(n) 1 \ n 

For the r.h.s. there is the well known formula [2], p. 27, = II ( 1 — — I = / fJ-{d) —;, if n 

n \ Pj J ^ d 

j=l ^ d\n 

has the distinct prime factors pj, j = 1,2, ...,M(n) = A001221(n) = \sqfkset{n)\ (see above for the 
notion 'squarefree kernel'). The Mobius function p entered here (see A008683). The PIE proof of eq. 
(30), given in the appendix A, uses this form of (p{n). If one uses also e*'^ = —1 the rewritten eq. (29) 
becomes 

n—l / , \ 

cy{n,-l) = {-!)— \[ 2 cos ( TT- J , n = 2, 3, .... (31) 

This can also be related to a special Sylvester sequence. In general the complex cyclotomic Sylvester- 
numbers Sy{a,b;n) are related to the sequence with three term recurrence /„ = a fn-i + b fn~.2, 
which has characteristic polynomial — ax — b with zeros a = a{a, b) = {a + Va2 + 4 6)/2 and 
p = j3{a,b) = (a - + Ab)/2. The definition is (see e.g., [19]) 

n-l 

Sy{a,b;n) := J] {a - (3 e^^'^) , n = 2,3, ... (32) 

gcd{k,rL) — l 
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For n = 1 one takes Sy{a,b;l) = q — /?. For our purpose the values are (a, 6) = (0,-1), i.e., 
(a, f3) = (z, —i), and 



n-l 



5y(0,-l;n) := (-l)A^ J] (1 + e^-*.) , n = 2, 3, (33) 



gcd{fc — 1 



with Sy{0, —1; 1) = 2i. This can be rewritten, using again eq. (30) and e*'^ = —1, as 



n-l y , N 

5y(0,-l;n) = (-1)'^(") JJ 2 cos ( vr- J , n = 2, 3, .... (34) 



gcd(k,n) — l 



Therefore 

5y(0,-l;n) := (-1)— cy(n, -1) ,n = 2, 3, .... (35) 

One could also include the n = 1 case, with ^{1) := 1 if one takes , . instead of the prefactor. 

After these preliminaries back to the number C(n, 0). 
Proposition 6: C(2m, 0), m G N 

C(2m,0) = "^^^""^ Syi0,-l;2m) = (-1)^ + ^ cy(2 m, -1) , mGN. (36) 

Proof: Use eqs. (19) and (35). The number of factors in eq. (19) with n = 2 m is , the degree 

5(2 m). Note that the restriction in the Sy{0, —l;2m) product is gcd{k,2m) = 1, while in the C(2m, 0) 
product it is gcd{k,Am) = 1, but this just says that only odd ks can contribute in both cases, and those 
odd numbers < 2m — 1 dividing 4m or 2m have to be omitted. Both restrictions exclude the same odd 
numbers. 

Corollary 1: C(2p, 0), p a prime 

p-i 

for prime p : C(2p,0) = (—1) 2 p. (37) 

Proof: 

This follows immediately from eqs. (19) and (31) if one uses the known formulae cy{2p,x) = cy{p, —x), 
cy{p, 1) = p, and evaluates the ip functions to obtain the sign. This formula follows from cy{n, x) = 
(x*^ — l)^^d) where the multiplicative Mobius function (with fi{p) = — 1) enters (see, e.g., [9], exercise 

d\n 

50 b., solution p. 506). 
Proposition 7: C(2'",x), m e No 

C{2"',x) = i{2'^-\x) , m e N, (38) 

with the integer polynomials t defined in eqs. (2) and (3) in terms of Chebyshev T- or S"— polynomials. 
Proof: Use eq. (16) and the Note added to the link [13], where it was shown in eq. (8) that 
22—1^1^2™+!, I) = 2r(2™-\|) =: t(2™-Sx), for m = 1,2,.... It is clear that C7(l, x) = x + 2. 
This proposition will later be generalized in theorem lA. 
Corollary 2: C(2™,0), m G No 

From the proposition 7 and the known fact that T(2n,0) = (—1)" follows that 

C(1,0) = +2, ^(2,0) = 0, C(22,0) = -2, C(2™,0) = +2, m > 3 . (39) 

Note that a standard formula for the cyclotomic polynomials cy{2"^,x) which involves divisors (see e.g., 
[9], p. 149, Exercise 50 b, with = cy{m,x)) leads to undetermined expressions for cy(2™,— 1) if 

m > 1. 
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For the following theorem lA on a formula for C{n,x) for even n, based on the divisor product represen- 
tation (dpr) of numbers [14], we need two lemmata. 

Lemma 3: Pairing in dpr(2 m) 

In the divisor product representation of an even number 2 m, m G N, for each o-factor in the numerator 
(resp. denominator) one finds exactly one a- factor in the denominator (resp. numerator) with argument 
ratio either 1 : 2 or 2 : 1. 

Proof: 

Due to the *-multiplication property of dprs (see the theorem in [14]) it is sufficient to prove this 
lemma for dpr{2pi ■■■ pjy) with odd primes pi, pN- Indeed, if the prime factorization of 2m is 
2^+1^61 with A; G No and odd primes pj and positive exponents e^, for j = 1, N, then 

dpr{2m) = (2^p^^~^ • • -p^"^) * dpr{2pi ■ ■ -pn) and the pairing property of the latter dpr will be pre- 
served after multiplication of each argument. Now the proof of the pairing property for dpr{2pi ■ ■ -pn) 
becomes elementary, once one splits the products in eq. (7) of [14] (note that N in this reference is now 
A^-|- 1) into those a-factors which are even {i.e., contain the even prime 2) and those which are odd (i.e., 
those composed of only odd primes). Indicate a-factors with odd arguments, being the product of k odd 
primes, by ao{.{k).). Then any a-factor in the numerator (resp. denominator) with some even argument 
2n is found in the denominator (resp. numerator) exactly once with the argument n. This is because for 

each factor in the numerator product na(2.(A^ — 2j).), with j £ {I,---, 



N 



}, say, a{2pi^ ■■■p, 



na(2.(iV-(2i + l)).), with j G {0,, 



*iV-2 j 1 1 

there is exactly one factor in the denominator product 11 ao(.(A^ — 2j).), namely the one which uses 
the same N — 2j odd primes a{pi-^ ■ ■ ■Pif^_2j)- Similarly, for each a-factor in the denominator product 

' N -I 

— ^ — }, there is exactly one factor in the numerator product 

Iiao{.{N — [2j + 1)).). The first factor in the numerator, a{2pi ■ ■ ■ p^) (j = 0) is paired with the single 
odd argument a-factor of the first product in the denominator Ilao{.{N).) = a{pi ■ ■ ■ pjy) ■ 
Note, that the example for the primorial A002110 (5) = 2310, given in [14], table 2 shows that the 
pairing occurs in general not between numerator and denominator factors at the same position, given the 
ordering prescription indicated by O (falling arguments) in eq. (7) in [14]. See the fifth position in the 
N = 5 case there. 

Recall also, for later purposes, that because dpr {2 m) has the same number of factors in the numerator 
and in the denominator, viz 2^ ( this balance holds true for all dpr{n), n > 2, due to proposition 5 in 
[14]), after the split into even and odd arguments the number of pairs with even argument in the numera- 
tor matches the one with even argument in the denominator. This will later lead to a balanced formula for 
C(2 m, x) in terms of f-polynomials (the same number of ts in the numerator and in the denominator, viz 
2^-1) . □ 

Lemma 4: 



t(f + l,x) - t(f - l,x) S(f + l,x) 



n x\ . /n 



2^ T' o = MT'^ , forn = 0(mod4) . (40) 



+ l,x) - t(f - l,x) 5(f - l,x) V4' 2J \4' 

This identity involves the integer polynomials t introduced earlier in eqs. (2) and (3). 
Proof: 

Use the known identity for Chebyshev polynomials, [20], p. 261, first line, specialized, T{m + l,x) — 
T{m — l,x) = 2 (x^ — l)U{m — 1,2;), written for x replaced by — , with S{m — 1,3;) = U{m — 1) ■^) 
and the definition of the t-polynomials from eq. (2). For n = Q{modA) this identity can be used 
in the numerator as well as in the denominator. The n— independent factor (^)^ ~ 1 drops out if 
x 7^ ±2. However, the lemma holds also for these x— values as can be seen after applying I'Hopital's 
rule, using the well known identity T'(x,n) = nU{n — l,x). In the second to last step the well known 
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identity, [20] p. 260, last line, written as 2r ^m, — ^ S{m — l,x) = 5(2 m — l,x) has been employed. 
□ 

Theorem lA: C(2m, x), m > 1, in terms of t-polynomials 

With the prime number factorization of 2 m = 2^ jy^ ■ ■ ■ with /c G N, odd primes "Pj and positive 
exponents ej, j = 1, A^, one has 

£(2^-1 v\^--- p%^,x) n £(2^-1 ,(iv - 2), x) n £(2^-1 ,(jv - 4)., x) • • • 

C 2m,a; = — - — ; — ; . 41 

^ ^ nt(2*^-i,(Ar-l).,x)nt(2'=-^(iV-3),x)--- ^ ' 

with i{2''-^ 40).,x) = t(2^~^Pr~V2'~^ • • •PaI^~\^)> and the products Ut{2''~'^ ^{K).,x) are over the 
("^) factors with K primes from the set {pi,P2, ■■■,Pn} multiplied by Pi^~^P2^~^ • • -p^"^, i-e., K- 
products of the form i{2^~^ p'!^~^ p'2~^ • • ■ p'^j^'^ Pii Pi2 " " 'Pik^^)- We used *(A;). instead of .{k). for the 
indices of the t polynomials to remind one of the extra factor (without the powers of 2 because the even 
prime has been extracted everywhere) due to the *-multiplication. 

Before we give the proof an example will illustrate this theorem. 
Example 3: k = 1, N = 2 

t{3-^ 5, X) t(3^ 5"^, X) t(135, x) t(225, x) 

Proof: 

This is based on [14] applied for the minimal polynomials "if (4 m, |) (see the definition for C(2m,x) 
given in eq. (16)). Let 2^^^ p^^ ■ ■ -p^^ be the prime number factorization for 4 m. Only the squarefree 
kernel 2pi- ■ -p^ is important due to the theorem in [14] for the divisor product representations (dprs) 
and the *-multiplication. This means that one has to multiply after the computation of 'I'(2pi • • -pn), 
using proposition 1 of [14], the first argument, the index in conventional notation, of each factor t{ni,x) 
in the the numerator, and of each t{mi,x) in the denominator (see [14], eq. (1)) with the number 
2^^61-1^62-1 _ _ _^ejy-i_ Here the pairing lemma 3 is crucial which carries over to the indices of the t- 
factors in the numerator and denominator of [14], eq. (1). Because there is always at least one factor of 
2 in the number which multiplies every t index after the *-multiplication, the pairing will occur always 
between indices which are 0(mod4) and 0(mod2). Then only the n even alternative in the definition 

of ^ from eq. (2) of [14], rewritten here as eq. (18), is relevant. Observe that the prefactors 

in this definition of t ^n, — ^ are not of interest, provided we take 2 ^— + l)'^^ ~ ^^(2^ ~ ^^2 
which is the monic integer polynomial i ^— + l,x^ — i ~ l)^;^. See eq. (18), and the discussion in 
connection with lemma 2. C(2m,x) is written here as a rational function of monic integer polynomials. 
For each replaced t ^n, — ^ j t ^— , — ^ (either in the numerator or denominator, depending on where the 

larger index appears), one can apply lemma 4- This is how for each such t-quotient one obtains i (^"^j ^ 

Here the fact that the larger index of a pair is always equivalent to 0(mo(i4) is important. From the 
structure of the numerator and denominator of the original dpr(2m) with the separation of the even 
and odd a-arguments (which carries over to the f-indices), discussed in the proof of lemma 3, one now 
finds the numerators and denominators of the theorem. Just search the products for the even i-indices. 
E.g., Ui{2''-\{N - 2). , x) in the numerator of the theorem originates from the quotient of products 
nt(2.(jV-2).,f) 



before the *-multiplication has to be applied. This leads, with eq. (18) and lemma 
/2fc+i \ 

4, toUi I ^ *(A^ — 2).,x j after multiplying each i index in the product with 2^ j5^^~^p2^~^ • • ■p'j^~ , 
where the * reminds one to multiply each index with this number divided by 2^^. 
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As announced at the end of the proof of lemma 3 the number of t polynomials in the numerator is 
the same as the one for the denominator, viz 2^"^, which is due to the sum in the Pascai-triangle A007318 
row No. N over even, resp. odd numbered positions. 

Note that this theorem lA, evaluated for x = leads in general to undetermined expressions, remember- 
ing that t{n, 0) = if n is odd, and 2 (— 1) 2" if n is even. However, a correct evaluation (using rHopitaPs 
rule) has to reproduce the result known from corollary 2. 

The following factorization of the monic integer i-polynomials is related to theorem lA. 
Theorem IB: Factorization of t-polynomials in terms of the minimal C-polynomials 

i{n,x) = Yl C{2n/d,x) = JJ C7(2'=+^ d, x) , (43) 

d\op{n) d\op{n) 

with op{n) = A000265(n), the odd part of n, and 2^ is the largest power of 2 dividing n. The exponents 
are k = k{n) = A007814 (n). /e G No . 

Before the proof we give an example. 
Example 4: n = 10 

-2 + 25x2 - 50 x^ + 35 x^ - 10 x*^ + x^° = t(10,x) = C(20, x) C(4, x) = 

(x^ - Sx*^ + 19 x^ - 12x2 + l)(x2 - 2) . (44) 

Proof: This is modeled after a similar proof in [30]. 

It is clear that both sides are monic integer polynomials. See the definition of t in eqs. (2), (3) and 
proposition 2. 

In order to check the degree we consider the even and odd n cases separately. If n = 2'^^"^ op{n) 
with k{n) = A007814(n) > 1 then ^ 5(2^'+^ d) = ^ known degree 5 of the 

d\op{n) d\op(n) 

C-polynomials, here for index > 1. Due to well known properties of the Euler totient function (see e.g., 

2fc{n)+2 _ 2fc(n)+l 

[2], Theorem 2.5, (a) and (c), p. 28 and Theorem 2.2., p. 26) this leads to op{n) = n, 

the degree of i{n,x). In the odd case, if n = op{n) then ^ ~ ~ ^S^i'^ degree of 

d\n d\n 

t{n, x). 

For the proof one compares the zeros of both sides. The zeros of iin,x) are x|"^ = 2 cos {(21 + 1) — ) 

V 2nJ 

for Z = 0, 1, n — 1. This is known from the zeros of the Chebyshev T-polynomials. The zeros 
of the (7(2'^"'"^ d, x)-polynomials are used in the form given by eq. (20): 2 cos ^(2 1' + 1) 2k+i ^ ) 

I' = 0,1,...,2'' d - 1, where gcd{2l' + 1,2^+^ d) = 1, i.e., gcd{2l' + l,d) = 1. Because the degrees 
match, it is sufficient to show that each zero of i{n,x) occurs on the r.h.s.. For n = 2^{n)op{n), with 
k{n) = A007814 (ra). consider Qcd(2l + l,2n) = gcd{2l + 1,2^+^ op{n)) = gcd{2l + l,op{n)) = g (some 

2 Z + 1 

odd number). IIence2/-|-l = (2 /'-|-1) (7, with some /' and op(n) = dg,i.e., d\op(n). Therefore = 

^ 2 re 

. Thus the t(re, x) zero 2 cos I (2 / -|- 1) - — I appears on the r.h.s. as one of the C{2^^^ d, x) zeros. 



2''+^d' ^ ' ' V ' 2n 

□. 

Remark 3: Derivation of Theorem lA from Theorem IB 

Theorem IB can be used as recurrence for the C-polynomials in terms of the t-polynomials. This is 

2 vr 

similar to the case treated in [30] for the minimal polynomials ^ of — (see A181875 / A181876 for their 
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coefficients). The solution of this recurrence has been given in [14]. Indeed, theorem lA has been derived 
above from this solution. Originally we found theorem IB starting from theorem lA building up iteratively 
a formula for t(2^~^p^^ • • -p^^ ,x) in terms of the C-polynomials. We give this formula as a corollary. 

Corollary 3: t-polynomials in terms of C-polynomials 

ei 62 ejv 

i{2^~^pl^...p'^,x) = n n n C'(2^pf •••Pr,^) ,A: = N, G No . (45) 

qi=0 92=0 ijjv=0 

In order to find a simplified expression for C(2 m + 1, x), m G Nq we need the following lemma in order 

t{even, |) 

to rewrite the quotient —. — , , ,,. . , given the definition eq. (18). 

t[odd, |) 

Lemma 5: 

for n = 2M + 1,M G Nq one has 

^(2n,f) ^ oM-n + 1^1)- ijr^ -!.§) _ 1 i^' - 4) S{n - 1, x) 
i(n,f) t(M + l,f) - £(M,|) 2*^+1 (x - 2)S(2MV2 + x) ^ ' 

X + 2 S{n-l,x) 

2M+1 5(n-l,V2"T^) ■ ^ ' 

Here new important monic integer polynomials enter the stage: 
Definition 1: q-polynomials 
With n G No define 

S(2n,x) ,, , 

6(2n, v2 + X) 

That this defines indeed monic integer polynomials of degree n is shown by the next lemma. 
Lemma 6: 

q(n,x) = (-l)"S(2n,V2 - x) = S{n,x) - 5(?i-l,x) . (49) 

Proof: 

It is known from the o.g.f. of the Chehyshev S-polynomials that the bisection yields S{2 n, y) = S{n, — 

2) + S{n — 1,7/^ — 2), or 5(2 n ,^/x + 2) = S{n,x) + S{n — l,x). With x replaced by —x one has 

5(2 n,\/2 — x) = (— l)"(5(n,x) — S{n — l,x)) which explains the second equation of the lemma. 

Therefore one has to prove S(2n,x) = S{n,x)^ — S{n — l,x)^. This identity can, for example, be 

proved using the o.g.f. for the square of the 5-polynomials (for their coefficient table see A181878 , 

also for the paper [15] given there as a link). This computation was based on the Binet-de Moivre 

I + z 1 

1 - Z 1 + (2 - x2)z + z2 

for {5(2n,x)}^g is (1 + z)/{l + (2 — x'^) z + z^) (from the bisection). Because the o.g.f. for 
S{n, x)^ — S{n — 1, x)2 is (1 + z) times the one for S{n, x)^ the completion of the proof is then obvious. 
□ 

Remark 4 : O.g.f. and coefficient array for the q-polynomials 



formula for the 5— polynomials. This o.g.f. was found to be tt- ^. The 0.17./. 



From lemma 6 the 0.17./. Q{z,x) := q{n,x)z^ = (1 — z)/{l — xz + z^) from the known o.g.f. 

n=0 

of the S'-polynomials. This shows that the q{n, x) coefficients constitute a iJiordan-array (infinite lower 
triangular ordinary convolution matrix) , which is in standard notation i , j , meaning 

i — X / X \ 

that the o.g.f. of the column No. m sequence is k ( k I . This is the triangle A130777 

1 + x'^ V 1 + / 
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where more information can be found. For example, the exphcit form for the coefficients is Q{n,m) = 

/n±m\ 

— 1) 2 ^ ifn>m>0 and otherwise. 

\ m J 

The pairing lemma 3 wih also be used in the proof of the following theorem. 
Theorem 2A: C(n, x), n > 3, odd, in terms of q-polynomials 

With the prime number factorization of n = • • • with odd primes pi, ■■■,Pn and positive exponents 
ej, j = 1, A^, one has 

n(llizl^\T\ n ( ^\ n „ ( '^dSPliIlzIlsIlilzl A . . . 

^, . 'i\ 2 ^-^J \.Ui<i2 'iy 2 '-^ ] Iln<i2<i3<i4 I 2 '-^i , , 

C[n,x) = J— I — j-^i ^- — . (50) 

iiii 2 '-^y iiii<j2<i3 y\ 2 ' 

Here each index ij of the products runs from 1 to N. 

Before we give the proof an examples will illustrates this theorem. 

Example 5: n = 3^ • 5 = 45 

This checks. 

Proof: This is analogous to the proof of theorem lA. Again the dpr(2n) representation, from which one 
derives ^{2n,x), with the prime number factorization for the odd n, is considered. The *-multiplication 
property allows to consider dpr[2pi - ■ ■ p]\j) with a subsequent multiplication of all arguments in the 
numerator and denominator by p'^~^ • • •p'j^^^- Because of the paired numerator/denominator structure 

t{2k,^) 

due to lemma 3 one finds, either in the numerator or in the denominator quotients of the type — — — -r- 

t(A:, f) 

fc + 1 / Tl — 1 \ 

which are up to a factor 2~2~ (x + 2) equal to q i — - — ,xj. Factors of powers of 2 are irrelevant 

(they have to cancel) because on both sides of theorem monic polynomials appear. The factors of x — 2 
also cancel because the number of g-polynomials in the numerator and denominator have also to match 
(see the remark at the end of the proof of lemma 3). This number is 2^"-"^. The structure of the 

numerator of the theorem originates from products with the even indexed t (^2k, — ^ polynomials in the 

numerator after *- multiplication. Before this multiplication on has in the numerator t ^2pi • • -pN, 

and t— products over all possibilities to leave out 2, 4, ... of the odd primes from the set {pi, pn}- The 

^-multiplication then leads to t (ji, — ^ and products over all possibilities to divide n by these 2,4, ... 

odd primes. Together with the pairing partners from the denominator this leads to the q-polynomials 
given in the numerator of the theorem. A similar argument produces the denominator g-polynomials. 
The number of these g— polynomials in the numerator, viz 2^~^ , matches the one for the denominator. 
□ 

Theorem 2B: Factorization of q-polynomials in terms of C-polynomials 

q{n,x) = Yl C{d,x) , n e N . (52) 

l<d|(2n+l) 



Example 6: n = 17 

This checks. 
Proof: 



/(17,x) = C{5,x)C{7,x)C{35,x) . (53) 
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This is again modeled after a similar proof in [30]. Compare this with the proof of theorem IB. 

It is clear that both sides are monic integer polynomials, and the degree fits due to n = ^2~^ ~ 

d\{2n+l) 

— (p{l)- See the properties of the Euler totient function mentioned above in the proof of theorem IB. 

The zeros of q{n,x) are = 2 cos — ^Y^' ' ~ 0,1,..., n — 1. This follows from definition 1 

and lemma 6 with the zeros of S{2 n, \/2 — x ) which are known from those of Chehyshev S'-polynomials. 
To show that each of these zeros appears on the r.h.s. for C{d,x) with some d\{2n + l),d 7^ 1, i.e., as 

2 cos f vr ^^-ii ) for some /' G {0, 1, ^^—^} and gcd(2l' + l,d) = 1 (see eq. 20), let gcd{2l + 1, 2n + 

\ Q J 2 
1) = c,, with some odd g with 2n + 1 = dg. Then 2/ + 1 = (2 /' + 1) 5 = (2 + 1)^^, and for each 

I G {0, 1, n - 1} there is one /' G {0, 1, ... , ^}. □ 
Remark 5: Derivation of Theorem 2A from Theorem 2B 

Theorem 2B can be used as recurrence for the C-polynomials in terms of the g-polynomials. See the 
remark 3. We give the solution in the following corollary. 

Corollary 4: q-polynomials in terms of C-polynomials 

With the prime number factorization of 2n + 1 = ■ ■ ■ , with odd primes, n G N, one has 

Qin,x) = n n C{pl'---P%n/Cil,x) . (54) 

qi=0 qN=0 

The division by C(l,x) = x + 2 was necessary because not all Qj-indices were originally allowed to 
vanish. 

These rational representations of C{n,x) do not lend itself to extraction of the value C(n,0) because 
undetermined — quotients appear. In the following we give the absolute term of C for prime indices. 

Propositions: C(p,0) 

For n = p, a prime, one has 



C{p,0) 



if p = 2 , 

(-1)^ if p = l(mod4) , (55) 
(-1)^ if p = 3(mod4) . 



Proof: In eq. (31) with n = p there is no gcd-restriction on the product. Therefore, one can use a 
known formula (see appendix B for a proof) 

n-i ^ (— 1)^~ if n is odd, 



n2cosfvr- = (56) 



k=i \ {) if n is even, 



to obtain cy{2, — 1) = (which is also clear from the definition: cy{2,x) = x + 1), and for odd primes 
cy{p, —1) = (—1)^"^ = +1. Because all k contribute in the product, one can use for even k the formula 

cos 1 vr ) = — cos ( TT ) which shows that one generates again all odd k contributions, how- 

\ P J \ P J 



p-3 

ever each with a minus sign. This leads to cy{p, — 1) = (—1)^"^ ( 2 cos ( vr 

1=0 ^ ^ 



p 
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Thus C{p,0) = ±1. The sign of C{p,0) is (— from eq. (19) for n = p. This number d+{p) of 
positive zeros has been found in proposition 5 to be if p = 1 {mod 4), and if p = 3 {mod 4). For 
C(2, 0) = 0, the sign is, of course, irrelevant. □ 

We close this section with a conjecture on the discriminant of the C— polynomials. The discriminant 
of a monic polynomial p of degree n can be written as the square of the determinant of an n x n Van- 
dermonde matrix Vn{x^^\ Xn"^) with elements {Vn)ij ■= {xf^^y , = 1, n and j = 0, n — 1 

with the zeros x^""^ of p . Here the 6{n) zeros of Cn are given in eq. (20). Another formula for the 
discriminant (see e.g., [26], Theorem 5.1, p. 218) is in terms of the derivative C'{n,x) and the ze- 

ros of C: (— 1) 2 I I C {n,x- ). The result is the sequence Discr{C{n,x)) = A193681 (ra) = 

i=l 

[1, 1, 1, 8, 5, 12, 49, 2048, 81, 2000, 14641, 2304, 371293, ...] for n > 1. The following conjecture is on the se- 

quence q{n) = — which is A215041 , [1, 2, 3, 2, 5, 3, 7, 2, 9, 5, 11, 9, 13, 7, 45, 2, 17, 27, 19, 25, ...1. 

Discr[C{n, x)) 

Conjecture: explicit form of the q-sequence 

0) q{l) = 1 (clear). 

1) If n = 2^ for /c G N then q{n) = 2. 

ii) \i n = p^ for odd prime p and A; € N then q{n) = p^^ 

iii) if n = 2'^'^ p{ii)^^i ■ ■ ■ p{iN)'^''^ with k2 € No, the ij — th odd primes p(ij), where 2 < ii < ^2 < 
... < IN, with TV G N if A;2 / and iV > 2 if /c2 = 0, then 

AT _ _ N 

q{n) = ^P(bO^'^"^^(^^)'""■■■^(^^^'"""^^"^'^■^with P(Ar,,) = \{{p{^l)-l), 
j=i 1=1,1^3 

= YY P"'^ ■ (57) 

odd p\n 

The last eq. follows from the degree S{n) = — ^. This last formula does, however, not work in the 

cases i) and ii) . One can compare this formula with the proven one for the discriminant of the cyclotomic 

polynomials (the minimal polynomials of exp{2Tr — ) (or any of the primitive n-th root of 1), as given in 

n 

[25], eq. (1) p. 297. For this (slightly rewritten) formula see also A004124 and A193679 . 
Example 7: 

ii) p = p(4) = 7,k = 3: q{7^) = g(343) = 7(^'+^)/2 ^ 1341068619663964900807. 

iii) n = 2^-3'^-7, q{n) = 32' -s-e 722.3 2 ^ 4316018525852839090954658176626149564980915348463203041. 
These values have been checked with the help of Maplel3 [18]. 

4 Splitting field Q{p{n)) for C(n,x), field extension and Galois group 

The algebraic number p{n) = 2 cos n £ N, with minimal polynomial C{n,x) over Q of degree 

S{n), has been studied in sect. 3. Each of these polynomials, being minimal, is irreducible. All roots 
have been given in eq. (19) (or eq. (20)). C{n,x) is also separable because all of its roots are distinct. 

Because 2 cos ( — ) = i{k,p{n)) (see eq. (2), and the coefficient array A127672), with the monic 
\ n J 

integer t-polynomials, each zero can be written as integer linear combination in the vector space basis 
< 1, p(n), p(n)^, p{nY^^'^~^ >, called the power basis. One has to reduce in t{k,p{n)) all powers p{n)P, 
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P > S{n) with the help of the equation C{n,p{n)) = 0. See Table 4 for the zeros of C{n,x), n = 1, ...,30 
written in this power basis. 

Example 8: n = 8 (octogon) 



The 6{8) = 4 zero s of C{8,x) = - 4x^ + 2 are, with p = p{8) = \/2 + y/2, ±/3 and ±(- 3p + 
p^) = ±\/ 2 — 1/2- In this case the degree 4 coincides with the number of DSRs in the upper half plane 
and the negative real axis. 

This shows that the extension of the rational field Q, called Q{p{n)), obtained by adjoining just one 
algebraic element (called a simple field extension) is the splitting field (Zerfallungskorper in German) 
for the polynomial C(n, x). Note that even though the polynomial C{n,x) is from the ring Z[x] one 
needs Q{p{n)) with rational coefficients rj for the general element a = Yl^j=o ^ rjp[ny. For example, 

p{8) = 2/9(8) — 2'°(^) ■ So™^ references for field extensions and Galois Theory are [6], [12], chpts. V 
and VI, [3], and the on-line lecture notes [16]. The dimension of of Q{p{n)) as a vector space over Q is 
6{n). This is the degree of the extension, denoted usually by [Q{p{n)) : Q], and it coincides with the 
degree of the minimal polynomial for p{n). Of course, it is a proper extension only if 5{n) > 2, i.e., for 
n > 4. This extension of Q is separable, i.e., the minimal polynomial for the general algebraic number 
a given above is separable (has only distinct zeros). It is a normal field extension, meaning that every 
irreducible rational polynomial with one root in Q{p{n)) splits completely over Q{p{n)). See, e.g., [6], 
Theorem 5.24, p. 108. 

We now consider a subgroup of the group of automorphisms of Q(p(n)), called Aut{Q{p{n)), which 
consists of those elements a which leave the subfield Q pointwise invariant (fixed point field Q): 
a : Q{p{n)) Q(p(n)), a ^ a (a), with a {/3) = /3 for ah /? € Q. The subgroup of ^Mt(Q(p(n)) of these 
so-called Q-automorphisms is called the Galois group of Q{p{n)) over Q, denoted by Qal{Q{p{n))/Q). 
Occasionally we abbreviate this with On- In order to find the elements a (we omit the label n) of this 
subgroup it is sufficient to know a{p{n)), because of the usual rules for automorphisms: (i) a{a + /5) = 
a{a) + cr(/3), (ii) a{al3) = a{a)a{P), and (Hi) a / =^ a{a) 7^ 0. Indeed, because of cj(l) = cr(l^) = 
(t(1)(t(1), one has a{l) = 1, and the images of products of p(n) are obtained from products of a{p{n)). 
Applying a on the equation C{n,p{n)) = (minimal polynomial), leads to C{n,a{p{n))) = 0, because 
the integer (rational) coefficients and are invariant under the Q- automorphism a we are looking for. 
Therefore we have exactly 6{n) distinct Q-automorphisms aj, j = 0, 6{n) — 1, determined from the 
distinct roots of C(n, a{p{n))), viz aj 

(pH) = with the zeros of C{n, x) ordered like in eq. (19) with 

increasing k values (see the Table 4). By the same token C{n,a{^'f'^)) = 0, for j = l,...,6{n) for every 
a. Because all roots are distinct (separable C) this leads to an isomorphism between Gal{Q{p{n))/Q) 
and a subgroup of the symmetric group Ss[n) ■ The Q-automorphisms a can therefore be identified with 

permutations of the roots of C (see e.g., [6] ch. 6.3, pp. 132 ff. One identifies the roots .fj"^ with j, and 

(jj with the permutation ( ^ ['^j(l) '^j(2) ••• '7j((5(n))]. This subgroup of S'5(„) is 

\aj[l) (Jj{2) ... aj[d[n)) J 

Abelian because only products of powers of p{n) appear and due to the automorphism property of a this 
carries over to the Galois group. 

Example 9: n = 5 (pentagon) Q-automorphisms 

a = ro 1 + ri p{5) with p{5) = ip, the golden section. In this case Q(v?) is as quadratic number 
field usually called Q(\/5) with the basis < l,(.p > for integers in Q(\/5) (see e.g., [10], ch. 14.3, p. 
207, where r = (93 — 1)). The two (5(5) = 2) Q-automorphisms are obtained from the solutions 

alip) = if ov = 1 — of C(5,n) = x^ — x — 1. Hence fio = id : cro(l) = 1, co(v) = V ^-nd 

ai : o"i(l) = 1, cri{ip) = 1 — (p. Because = ctq, ^5 is generated by cji, hence the Galois group 
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^aZ(Q((/j)/Q) is the cychc group (also known as additive group Z/2Z or C2, but we reserve C for the 
minimal polynomials). The fixed field for is Q; for the trivial subgroup with element it is 



Example 10: n = 7 (heptagon) Q- automorphisms 

The three zeros (5(7) = 3) of (7(7, n) are = />(7), = -1 - p(7) + p{lf and = 2 - p(7)2 
(see Tahle 4)- In the sequel we omit the argument 7. Computing powers of p modulo C(7, p) = one 
finds the Q-automorphisnis ao = id, ai{p = ^i) = ^2, o-i(l2) = Cs, (^liiz) = ?i and a2{p = Ci) = 6, 
<72('?2) = '^i, cy2{£,z) = '?2- Therefore, the identification with 53 permutations is (Tq = [12 3] = e, 
(Ji = [2 31] and (T2 = [312]. Each aj permutation can be depicted in a circle diagram with vertices 
labeled 1, 2 and 3 and directed edges, also allowing for loops. The only (Abelian) group of order three is 
the cyclic Z^, subgroup of ^3. 

We can also characterize these Galois groups Gn by giving their cycle structure and depict them as cycle 
graphs. For cycle graphs see, e.g., [31] "Cycle graph" and "List of small groups" with all cycle graphs for 
groups (also non- Abelian ones) of order 1, ... , 16. In order to manage powers of Galois group elements 
we first need some fundamental identities of Chebyshev T-polynomials, whence t-polynomials. 

We will need iterations of t polynomials governed by the following well known identity. 

Lemma 7: Iteration of t-polynomials 

t{n,t{m, x)) = t{nm,x), n,m E Nq . (58) 

Proof: See [26], Exercise 1.1.6, p. 5, first with the trigonometric definition of Ciiebysiiev T-polynomials 
which then carries over to the general polynomials defined by their recurrence relation. This identity is 
then rewritten for t-polynomials. □ 

Lemma 8: mod n reduction of t-polynomials in the variable p(n) 

i{k,p{n)) = (-l)LnJ £(/c(modn),p(n)), n G N, A; e Z . (59) 

(-l)L =: Pn{k), the parity of |_|J , will become important in the following. Sometimes |_^J is called 
quotient and denoted also by k\n. 

Proof: This follows, with k = In + r, trivially from the trigonometric identity cos ^— (In + r)^ = 

cos(^ttI + — ^ = (-1)' cos (^-r^, with r G {0, 1, ...,n-l}. □ 

To simplify notation we will use also tn{x) for t{n, x). From lemmata 7 and 8 we have e.g., t3{t3{p{7)) = 
ig{p{7)) = -i-2{p{7)) = -{p{7f - 2). This is also +h{p{7)) from the identity cos Q(n ± /)) = 

— cos (^—l^. This proves the following lemma. 

Lemma 9: Symmetry relation of t polynomials 

i{n-l,p{n)) = -i{l,p{n)), ?i G N, / G {0,1,..., n} . (60) 

Lemma 7 which used the trigonometric p{n) definition can be rewritten as a congruence for the t- 
polynomials with indeterminate x. This is because all what is needed is that p{n) is a zero of C{n,x). 

Corollary 5: Congruence for t polynomials modulo C-polynomials 

i{k,x) = (-l)LnJ i{k{modn),x){modC{n,x)) , n G N, /c G Z . (61) 

An example will illustrate this before we give an another proof of this corollary based on known T- 
polynomial identities. 

Example 11: Congruence for n = 7, k = 9 

i{9,x) = x^ - dx"^ + 27 x^ - 30 x^ + 9x and C{7,x) = x^ - x^ - 2x + 1. Polynomial 
division shows that t{9,x) = {x^ + x^ — 6a;^ — 5x^ + 9x^ + 5x — 2)C{7,x) + (— -|- 2), hence 
i{9,x) = -i{2,x){modC{7,x)). 
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The alternative proof of corollary 5 is based on the following factorization of Chebyshev S'-polynomials 
in terms of the minimal polynomials \I'(n, x) of cos ( — ) over Q (for these polynomials see A049310 



n 



and A181875/A181876 ). 

Proposition 9: Factorization of S-polynomials in terms of 'J'-polynomials 

5(n-l,x) = 2"-^ Yl *((i,|),nGN. (62) 

2<d I 2n 

As usual, the undefined (empty) product is defined to be 1. 

Proof: Start with the eq. (3) of [30], p. 471, written for n — > 2n and for the i-polynomials instead of 
the T-polynomials. This is t(n + 1, x) — i{n—l,x) = 2"^"*^ ^(^d,^^. Then use the identity [20], p. 

d\2n 

261, first line, with m ^ n, n ^ m and [/-polynomials replaced by 5— polynomials. This leads to the 
identity 

i{n + m, x) — i{n — m,x) = (x^ — 4) S{n — l,x) S{m — 1, x) , m < n S N . (63) 

Here we only need the case m = 1. Then dividing out 2^(1, |) = x — 2 and 2^'(2, |) = x + 2 leads 
to the claimed identity. □ 

To end the preparation for an independent proof of corollary 5 we state: 

Corollary 6: C-polynomial divides some family of S-polynomials 

C(n,x) |5(/n- l,x), n > 2, / G N . (64) 

This is clear from the definition of C in eq. (16) and the fact that 2 < 2n\2ln, for / > 1. E.g., 
C(10,x) = x^ - 5x2 + 5 divides the family {5(9, x), 5(19, x), 5(29, x), ...} 

Proof of corollary 5: From the corollary 6 and the identity eq. (63) with n ^ In and m ^ r one sees 
thatt(A:,x) = t{ln + r, x) = t{ln — r,x) {modC{n,x)). Nowi(/n— r, x) =t((/ — l)n+ (n — r),x), 
and one can use h -.= 1 — 1 and ri := n — r (remember that r € {0, 1, n — 1}, hence ri has identical 
range) as new I and r variables in this congruence, to get t{li n + ri, x) = t{li n — ri, x) {modC{n, x)) = 
t{{l — 2) n + r, x) {modC{n, x)). This can be continued until one ends up with 

t{On+ k{modn),x) {modC{n,x)). □ 

Now we are in a position to compute powers of elements of the Galois group Gn- The subset of odd 
numbers 2 / + 1 < n entering the product in eq. (20) will be denoted by Ai{n). There are 6{n) (degree 
of C(n, x)) such odd numbers. 

Definition 2: The fundamental set A4{n) 



M{n) := <^ 2 / + 1 



/ G <^ 
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andfi(C(i(2^ + l,n) = 1| = {mi(n), m5(„) (n)}, (65) 

with mi(n) = 1 and we use the order mi{n) < mj{n) if i < j. For n = 1 one takes Al(l) = {1}. 
Example 12: \M{2)\ = 6{2) = I, M{2) = {1}; \M{U)\ = 5{U) = Q , M{U) = {1,3,5,9,11,13}. 
For Ai{n) see the row No. n of the array A216319 . 

For later purpose we define here, for odd n, the extended fundamental set M(n) and its first difference 
set A M{n) 

Definition 3: The fundamental extended set 7W(n) for odd n 
For n odd, > 1: M.{n) := {0, mi(n) = 1, (n) = n — 2,n + 2} . 

Thus |A^(n)| = 6{n) + 2. Note that gcd{n ± 2,n) = 1 for odd n (proof by assuming the contrary: 
gcd{n + 2, n) = d > 2 because n and n + 2 are odd. Then ti | (n + 2) and d \ n, hence d \ ((n + 2) — n), 
d I 2, implying d = 2 or d = 1, but d > 2.) The reason for defining this extended set is that for odd n the 
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first difference set AA4{n) = {1, A?n2(n), Am^l-^) (n), 4} with Amj{n) := mj{n) — mj-i{n), will 
become important later on. Ami(n) = 1 and Ams(^n)+i{^) = 4 for each odd n. 
With definition 2 eq. (20) implies, 

CTjipin)) = crj{i["''^) = = 2 cos (^^?Tij+i(n)) = 

t(m,+i(n), p(n)) for j e {0, 1, 6{n) - 1} . (66) 

E.g., n = 7: 6{7) = 3, M{7) = {1,3,5}, ao(p(7)) = p(7), ai(p(7)) = t(3,p(7)), and a2(p(7)) = 
t(5,p(7)). 

Because t is a rational integer polynomial 

c^jipin)) = aj{aj{p{n))) = i{mj+i{n),aj{pin))) = £(mj+i(n), £(mj+i, p(n))) , (67) 

and with lemma 6 this becomes aj{p{n)) = i{{mj^i(n))'^ , p(n)). In general we have 

a^{p{n)) = ^((mJ^.l(n))^/5(n)), fori G {0,1, ...,6{n) - 1} ,n,A: G N . (68) 

Instead of powers of aj{p{n)) we can therefore consider powers of mj_|_i(n). Lemmata 8 and 9 are now 
employed to prove that a t polynomial with a product of elements from 7W(n) as its first argument (or 
index) is again a t polynomial with first argument from M{n). In this way one can build sequences of 
powers, starting from any element of Ai{n). Trivially, 1^ = 1. Before proving this closure of M{n) 
under powers, provided the rules for t polynomials are taken into account, we give two examples, and 
then define a new equivalence relation on the integers, called {Moddn), denoted by S . 

Example 13: Cycle structure for n = 12 (dodecagon) 

(5(12) = 4, A^(12) = {1,5,7,11}. 52 = 25 = 1 (mod 12), reflecting lemma 8. The sign p„(25) in eq. 
(59) is here +. (If later the sign will be — , the mod n result will be underlined.) The first 2-cycle 
is therefore [5,1]. Similarly, 7^ = 49 = 1 (mod 12) (sign +), whence the second 2-cycle is [7,1], and 
finally, 11^ = 121 = 1 [mod 12) (sign +), producing the third 2-cycle [11, 1]. This result appears as the 
n = 12 entry in Table 6. In this example the Galois group is not generated by one element, hence it is 
non-cyclic. In fact, Gal{Q{p{12))/Q) = Z2 x Z2 = Z2 (see second entry in Table 1). The corresponding 
cycle graph is shown in Figure 4 as the first entry, where the shaded (colored) vertex stands for 1 and 
the open vertices should here be labeled with 5, 7 and 11. The cycle structure is 23 (three 2-cycles). In 
this example it was not necessary to employ lemma 9 because the signs were always -|-. 
Example 14: Cycle structure for n = 7 (heptagon) 

5(7) = 3, M{7) = {1,3,5}. 3^ = 2 {mod 7), where now the sign — in eq. (59) is remembered by 
the underlining. 2 ^ -A4(7), and now lemma 9 is used to rewrite this 2 as 7 — 2 = 5. Therefore, 
3^ Z 5 G M.{7) (or 3^ = 5 {Modd7)). The symbol Z (or Modd 7) is used for the congruence in the new 
sense, due to lemmata 8 and 9. Then 3-5 = 15 = 1 {mod n) (sign -|-). The first 3-cycle is therefore 
[3, 5, 1]. Here 3 generates all the elements of (7), and Gal{Q{p{7))/Q) = Z3, the cyclic group of order 
3 = 6{7). The corresponding cycle graph is a circle with three vertices, one of them, labeled 1, is shaded 
(colored) and the other two open ones are labeled by 3 and 5. 

This brings us to the definition of an equivalence relation S (or Moddn) over the integers Z. Remember 
that the floor function for negative arguments is deflned as [—x\ = — [x\ if x G No, and [—x\ = 

-{[x\ + 1) if < X ^ n. 

Definition 4: Equivalence relation on Z 

For k,l G Z, n G N: /c S / CLn{k) = an{l), with the map an '■ ^ In := {0, 1, n — 1}, k 1— )■ an{k), 
where 

{rn{k) if pn{k) = +1 , 
(69) 
rn{-k) if pn{k) = -1 , 
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where we used the division algorithm to write k = qn{k)n + rn{k), with the quotient qnik) G Z and 
the residue r.„(/c) E /„. Note that g„(A;) = — . The sign p.„(/c) = ( — 1)L"J = (—1)'^"^'^^ corresponding 
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to the parity of Qnik), appeared already in lemma 8. a„(l) = because —1 = [modi). 

Instead of kJbl we also write k = l{Moddn) (this should not to be confused with modn). Therefore 
the sequence could also be called Modd n. The first of these 2n-periodic sequences are found in 
A000007 (ra + n.n > 0, (the 0-sequence), A000035 . A193680 . A193682 , A203571 , A203572 and A204453 , 
for n = 1, 7, respectively. 

The smallest non-negative residue system mod n, viz 0, 1, ...,n — 1 is used here. For the residue of k 
modulo n Maple [18] uses r„(A;) = modp{k,n). We also use k{modn) for rn{k). The reader should 
verify that Jh is indeed an equivalence relation satisfying refiexivity, symmetry and transitivity. The 
disjoint and exhaustive equivalence classes are given by {[0], [1], .., [n — 1]}, called the smallest non- 
negative complete representative classes (or residue classes) Moddn (we omit the index n at the classes: 
[m] = nlTn], written this way in order to distinguish this class from the ordinary one [m]n used in the 
arithmetic modn). These classes are defined by [m] = {/ G Z|/ Zim} = {/ G Z|a„(Z) = m}. Because 
fn{—k) = if /c = [mod n) and rn{—k) = n — rn{k) if A; ^ [mod n) (later listed as lemma If) one 
can characterize these residue classes also in the following way. 

Lemma 10: Complete residue classes Moddn 

/ G [0] <^ / = 0{modn). Equivalently, [0] = [0]„ , 

I G [1] ^ I = l{mod2n) or = -l{mod2n). Equivalently, [1] = [l]2n U [2n- l]2n , 

I G [n — l]<^l = n — 1 {mod2n) or = — (n — 1) {mod2n). 

Equivalently, [n-\] = [n - l]2n U [n + 1]2„ . (70) 

For example, take n = 7, then [3] = 7[3] = [3]i4 U [14-3]i4 = {..., -25, -11, 3, 17, ...} U 
{...-17,-3,11,25,...} = {...,-25,-17,-11,-3,3,11,17,25,...}. If n = 2 one has 2 [0] = [0)2 and 
2(1] = [1)2- The first differences in the class [0] are, of course, n, and in the class [m], for m = l,2,...,n — 1, 
they alternate between Ai = 2 (n — m) and A2 = 2 m, e.^., n = 7, m = 3, Ai = 8, A2 = 6 : 
3, 11, 17, 25, 31, .... This difference alternation invalidates certain theorems known for modn. E.g., 
Theorem 53 of reference [10], p. 50, is no longer true: take m = 5, n = 7, a = 3, 6 = 17. a and h 
belong both to the class 5(3] as well as 7(3] but they do obviously not both belong to the class 35 [3] . The 
following lemma 11 shows that is is sufficient to know the positive values, and append the negative of 
these values for each class [?n], for m G {0,l,...,n — 1}. 

Lemma 11: Antisymmetry of the classes [m], for m G {0, ...,n — 1} 

For n G N and every m G {0, ...,n — 1} the elements of the equivalence class [m] are antisymmetric 
around 0. 

Proof: This is obvious for the class [0] = {... — 2 n, — n, 0, n, 2 n, ...}. For m > the negative of every posi- 
tive numbers of m {mod2n) appears as a negative number of (2n— m) (mod2n), and vice versa: the nega- 
tive of every positive numbers of [2 n—m]2n appears as a negative one of [?7T.]2n- This is true because — (m+ 
I2n) = {2n — m) — (/ + 1) 2 n, for every I G N, and similarly, — ((2 n — m) + I2n) = m — (/ + 1) 2 n, for 
every / G N. □ 

This leads immediately to the following corollary. 

Corollary 7: Non-negative elements of the classes [m] 

For n G N and every m G {0, n — 1} one has 

[m]> = [m]2n,> U [2n-m]2„,>, and [m] = [m]> U -([m]>) . (71) 
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Here we used the notations [m]> and [•]2n,> or [•]2n,> to denote the subset of non-negative numbers 
of [m] and the non-negative or positive numbers of the ordinary residue classes mod2n, respectively. Of 
course appears only in the class [0]. — ([?7i]>) is obtained from the set [m] > (excluding in the case of 
class [0]) by taking all elements negative. Note that [—m] is not used here. 

Withg'(n,m) := 5(C(i(m, 2 Ji-m) one has, for jn G {1, 2, n-1}, [m]2„,> = g{n,m)[m/ g{n,m)]2n/g(n,m),> 
and [2n — m]2n,> = ^(n, m)[(2n — 'm)/g{n,m)]2n/g{n,m),>- This is. of course, only of interest if 
g{n,m) / 1. E.g., n = 6, m = 3, g{n,m) = 3: [3]i2,> =3* [1]4,> and [9]i2,> = 3* [3]4,>, where 
again k * \p]q,> is the set with all members of the set \p]q,> multiplied with k. This is obvious. 

The trivial formula for the members of the residue classes [m]>, considered as sequences of increasing 
numbers called {c{n,m; k)}^^, is 

{{k — l)n if m = , 

(72) 
[|j2n+(-l)'=m if m G {l,2,...,n- 1} . 

The nonnegative members of the complete residue classes {Moddn) for n = 3,4,5,6, and 7 are found in 
A088520 , A203575 , A090298 , A092260 , and A113807 . Sometimes has to be added, in order to obtain 
the class [0]. Of course, these complete residue classes can be recorded as a permutation sequence of the 
non- negative integers. 

We now list several lemmata (the trivial lemma 17 has already been used) in order to prepare for the 
proof of the multiplicative structure of these equivalence classes. 

Lemma 12: Parity of Moddn residue classes 

For even n the parity of the members of the residue class ri\jn\-, m G {0,l,...,n — 1} coincides with the 
one of m. If n is odd this is also true for the classes with m G {1, 2, ... , n — 1}, and for m = the parity 
of the elements alternates, starting with + (for even). 

Proof: The case of the residue class [0] is clear for even or odd n because of its members Omodn (see 
lemma 10). Similarly, for the other m values, because then mod2n applies. 

Lemma 13: Periodicity of the parity sequence pn with period length 2 n 

Pn{k) = Pn{k + 2n/) for / G Z, i.e., Pn{k) = Pn(?'2n(^)) = Pn{k {mod2n)), for /c G Z,n G N. (73) 

Proof: This is obvious from the definition oi Pn{k) given in lemma 8, eq. (59), with definition 4, eq. (69). 

For the second part use k = 2nq2n {k) + r2n{k). 

Lemma 14: (A)symmetry of sequence pn around k = 

{+Pn{k) \i k = 0{modn) 
(74) 
—Pn{k) , if A: ^ {){modn) . 

Proof: This follows immediately from the property of the \_—x\ function mentioned above before defini- 
tion 4- n 

Lemma 15: Product formula for the residue rn 

For n G N and k,l £ 7j one has: rn{kl) = rn{k) rn{l) (modn). 
Proof: Just multiply A; = qn{k) n + rn{k) with. I = qn{l) n-\- rn{l). 
Lemma 16: Residue r2n from rn 

r rnik) iSr2n{k) G {0, 1, n-1}. 

For n G N and A; G Z one has: r2n{k) = \ 

[ rn{k) +n iff r2„(A;) G {n, n + 1, 2n - 1}. 

Proof: Obvious for r2 n {k) G {0, ...,n-l} as weh as G {n, 2n-l}. □ 
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Lemma 17: Residue for negative numbers 

{0 if A; = {modn) , 

n — rn{k) if k ^ 0{modn) . 
Proof: Obvious. □ 
Lemma 18: Symmetry of the an (or Moddn) sequence around k = 

For n G N and A; G Nq one has: o„(— A;) = a„(/c). 

Proof: This is proved for the two cases k = 0{modn) and k ^ 0{modn) separately. In the first case 
fnik) = 0, as well as rn{—k) = , hence a„(A; = (modn)) = 0, which is symmetric around A; = 0. In 
the other case, we employ lemma 13, noting that r2n{k) 7^ 0, n (otherwise A; = {modn)). The two cases 
r2n{k) G {1, 2, n - 1} and r2n(^) S {n + 1, n + 2, 2n - 1} have p„(A;) = Pn{r2n{k)) equal +1 or 
— 1, respectively. Then with the second alternative of lemma I4 one finds an{—k) = rn{—h) ii pn{—k) = 
—Pn{k) = +1 and an{—k) = rn{k) if A;) = —pn{k) = —1 which coincides with the definition of 
an{+k). □ 

Now we turn to the arithmetic structure of the Moddn residue classes. It is clear from the following 
counter-example that addition cannot be done class-wise. Consider n = 6, A: = 2 and / = 7. Then 
a6(2 + 7) = 3 but a6{ae{2) + a^iJ)) = 05(2 5) = 05(7) = 5. In other words, 2 6,2 and 7^5, but 
2 + 7 = 9^3 but (2 5) = 7^5, and 3 is not equivalent to 5{Modd6). Similarly from (a - l)XiO 
does in general not follow a S 1, because for n > 2 in the latter case a can also be of the form —l + k2n 
if it belongs to a class [m] with positive m, whereas in the first case it has to be of the form 1 + k' n, 
which can only match for n = 1 and 2. However, it turns out that multiplication can be done class- wise. 
This is the content of the following proposition. 
Proposition 10: Modd n residue classes are multiplicative 
For n G N and k,l £ 1^ one has: 

an{kl) = an{an{k) an{l)) i.e., kl Jh an{k) an{l) , i.e.. A;/ = an{k) an{l) {Modd n) . (75) 

Before giving the proof, consider the example n = Q, k = 2 and I = 7. Now 05(2 • 7) = 2 and 
06(06(2) • 06(7)) = 06(2 • 5) = 05(10) = 2. Or stated equivalently 2 • 7 = 14 6 2 6 2 • 5 ^ 05(2) 05(7). 
Proof: 

i) Due to the symmetry of o^ (see lemma 18) it is clear that it is sufficient to consider only non- negative 
k and /. 

ii) Consider first the cases k = 0{modn) or I = 0{modn). If m = 0{modn) then an{m) = 0. This 
follows for both alternatives in eq. (69). Therefore, if A: = [modn), an{kl) = because kl = (modn) 
for every and o„(0 • an{l)) = o„(0) = 0, proving the assertion. In the other case, I = 0{modn), the 
proof is done analogously. 

iii) Now r„(A:) and r„(/) are non-vanishing, and k and I are positive. Four cases are distinguished 
according to the signs of {pn{k) , Pn{l)) , viz (+,+), (-,-), (+,-) and (-,+). 

(-|-,-|-): In this case, due to the 2 n-periodicity of Pn (see lemma 13), r2n{k) and r2n{l) are both from 
{1, 2, n — 1}, hence 

(+,+) : r2n{k) = rn{k) and r2n{l) = rn{l) ■ (76) 

Also, from eq. (69), one obtains in this case an{k) an{l) = rn{k) rn{l). There is an alternative for 
<in{Tn{k) rn{l)), depending on Pn{rn{k) rn{l)) being +1 or —1. Both cases are possible as the following 
examples for n = 6 show: pq{3-A) = pq{12 {modl2)) = pq{0) = +1 andp6(3-7) = pq {21 (mod 12)) = 
Pq{9) = —1. First the argument of Pn is rewritten with the help of eq. (76). Pn{i'n{k) rn{l)) = 
Pn{'f2n{k)r2n{l))- Due to the (2 n)-periodicity (ZeTOTOa 15) this is (r2 n (fc) ^2 n(0 "^od(2 n)) = Pn{i'2n{kl)) 
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p{k I), due to lemma i5 with n — )• 2 n, and again the (2 n)-periodicity. In the first alternative, the +1 case, 
O'Tiio-nik) an{l)) = rn{k)rn{l) {modn) = rn{kl), again from lemma 15. This is just a{kl) li pn{kl) = 
+ 1, proving the assertion for this alternative. In the other case, Pn{kl) = —1, an{'rn{k)rn{l)) = 
— {Tn{k)rn{l)){'modn) which is rewritten with mod n-arithmetic and lemma 15 as —rn{kl){modn) = 
~{kl){modn) = rn{—kl). This coincides with a{kl) for this alternative, proving the assertion. 

(-,-): Now we have from lemma 16 



Here Pn{{fn{—k)rn{—l)) is rewritten with lemma ^7, eq. (77), the 2n periodicity, and lemma 15, 
as follows. p„((n - r„(A;)) (n - r„(/))) = p„,((2 n - r2„(A:)) (2 n - r2„(/))) = Pn{r2n{k)r2n{l)) = 
Pn{r2 n(''2 n(^) ^2 n(0)) — Pnii"2 n{k) 1^2 nil-) {^od2 Ti)) . With lemma 15 (with n >■ 2n this becomes 
Pn{T2n{kl)) = Pnikl), again from the (2 n)-periodicity. In the first alternative Pn{kl) = +1 and 
an{fn{—k)rn{—l)) = rn{rn{—k)rn{—l))- With iemma i7and mod n-arithmetic this is r„((—r„(A;)) (—r„(Z))) = 
fn{fn{k)rn{l)), and with lemma ^ 5 this becomes rn{kl), coinciding with an{kl) for this alternative. For 
the other alternative, p„(fc = — 1, a„(r„(— /c) r„(— /)) = —{rn{—k)rn{—l)){modn). This becomes, with 
lemma 17, mod n- arithmetic and lemma 15 

~{Tn{k) rn{l)) {modn) = —rn{kl){modn). This vanishes if r„( A; /) = which means fcZ = 0{modn) 
(which is possible, e.g., n = 6,k = 2,1 = 3), and then this coincides with the claim which is for this 
alternative a{kl) = rn{—{kl)) = 0. If r„(A;Z) ^ then —rn{kl) {mod n) — n — T,fi{kl) which also 
coincides with the claim a{kl) = n — rn{kl) \ikl ^ 0{modn). 

We skip the proofs of the other two cases, (+, — ) and (— , +), which run along the same line. Here one 
arrives first at Pn{—{k I)), and in order to compare it with p{k I) both alternatives in lemma ^ 7have to be 
considered like in the just considered second alternative. □ 

For the computation of the cycle structure of the Galois group Qn = G<il{Q{p{n))/Q) we are only 
interested, due to eq. (20), in odd numbers relatively prime to n. Contrary to ordinary mod n-arithmetic 
where the set of odd numbers O := {21 + 1 1 / G Z} is in general not closed under multiplication 
{e.g., 5 • 5 = 25 = 4 (mod 7)), it will be shown that O is closed under Moddn multiplication. Of 
interest are the units, the elements which have inverses, in order to see the expected group structure. 
First consider the reduced set O*, given by the odd numbers relatively prime to n. The negative odd 
numbers in this set are just the negative of the positive odd numbers, therefore it will suffice to consider 
0;^> := {2^ + 1 1 / e No, 5'cd(2/ + l,n) = 1}. If n is a power of 2 the set wih be O, with the o.g.f. 

X 

G{x) = — j2 (1 + ^) fo^ the sequence of positive odd numbers {o2^y{k) := 2k + Ij^Lo- For the 

other even numbers n only the odd numbers relatively prime to the squarefree kernel of n, called sqfk{n), 
(see A007947 , encountered already several times), will enter the discussion. Therefore, besides the just 
considered (trivial) case of the even prime 2, the set O* is only relevant for squarefree odd moduli, either 
prime or composite. We consider first the case of odd primes n = p, and give the o.g.f. of the sequences 
of numbers from Op^>, called Opy{k), as well as an explicit formula in terms of floor functions. These 
are the positive odd integers without odd multiples of the odd p. 

Proposition 11: Odd prime moduli, o.g.f. and explicit formula for O* > elements 



r2n{k) 



n + rn{k) and r2„(/) 



n + r„(/) . 



(77) 



oo 



With Gp{x) := o*p^^{k)x^, for odd primes p, one has 



P-3 

2 



Gp{x) 



X 



< 1 + 2 ^ a;'' (1 + xV") + 4^ + x' 



.p-i 



(78) 



(1 - xP-i) (1 - x) 



k=l 



29 



and 



p,> 



(k) 



2k 



1 + 2 



if k = , 



if A; G N . 



(79) 



We have used the even number for k = such that the Gp{x) sum can start with A: = 0. 

Proof: This proposition will become a corollary to the later treated general case of odd squarefree moduli 
n in proposition 13. 

Example 15: O.g.f.s and formula for reduced odd numbers for modulus p = 7. 



Grix) 
07>(n) 



(l-x6)(l-x) 



{l + 2(x + x^) + 4x^ + 2(x^ + x^) + x^} 



2n 



1 + 2 



n 



6 



n > 1. 



(80) 
(81) 



The instances for p = 3,5,7,11,13 and 17 are found under A007310 , A045572 , A162699 , A204454 , 
A204457 , and A204458 , respectively. 

In order to prepare for the general case of odd squarefree modulus n, we state a proposition on the 
structure of the reduced odd numbers set O* 

Proposition 12: Mirror symmetry and modular periodicity of O* > for odd n 

i) Mirror symmetry. For k G {1, 2, 5(n)} one has: 



{2 6{n)-{k-l)) 



2n 



where o^^{k) = mk{n) from 7W(n), given in definition 2, eq. (65). 
ii) mod 2 n periodicity: For k £N one has: 



(k) 



{k + 2 5{n)) {mod2n) 



(82) 



(83) 



Written as a relation between neighboring fundamental units, numbered by > 1, this becomes the 
following statement. For k G {2 (A^ - 1) 5{n) + 1, 2N 5{n)], with A^ G {2, 3, 4, ...}, one has: 



{k) = {N -l)2n+ ol^{k - [N -l)25{n)) 



(84) 



5{n) is the degree of the minimal polynomial C(n, x) for the algebraic number p{n) introduced in section 2. 

u>(n) u){n) 

Note that if n = JJ^ pj with distinct odd primes p, andw(n) = A001221(n), then 2 5(n) = JJ(pj-l). 

i=i i=i 
L{n) := 2 5{n) is the length of the fundamental A^-units. 

Before we give the proof consider figure 3 for the case n = 3 • 5 = 15 with 5{n) = 4. The second 
statement ii) concerns the relation of the numbers of the second fundamental unit (A^ = 2) to the one 
in the A^ = 1 unit. E.g., the odd number 37 for A; = 2-4 + 2 = 10 is equal to 30 + 0^5 >(10 - 2-4) = 
30 + 0^5^(2) = 30 + 7, which checks. The statement i) shows the mirror symmetry within the first 
(and any other) unit of length L{l'5i) = 2 • 4 = 8. In figure 3 this symmetry is indicated by the brackets 
below the first unit, and it is a symmetry around the missing number n = 15. Missing numbers have 
been indicated by a dot. It is the pattern of missing odd numbers which is mirror symmetric, not the 
one of the actual values of the odd numbers. But the relation between the odd numbers in the second 
half of a unit and the first one then follows, and is given by the statement of the proposition. In figure 
3 P{n) := 2n = 30 is the shift for the o> „ values from the A^ = 1 to the N = 2 unit (or any of the 
neighboring units), and p{n) := n + 1 = 16 is the shift for these values from the first half of every unit 
to the second half. 
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Figure 3: Structure of the sequence of mod n reduced odd numbers, n=15 
N=l N=2 
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Proof: i) It is clear from the degree of the minimal polynomial and eq. (20) that the number of modn 
reduced positive odd numbers smaller than 2n is 5{2n) = 2 6{n) (from the definition of 5{n) in terms of 
Euier's (p function). We now determine 5{n) reduced odd numbers modn which lie between n and 2n, 
by mirroring the 6{n) elements of M{n) around the position where the missing number n is situated, 
which is between the position k = 6{n) and the next one. The mirror symmetry refers to the gaps. The 
number mk{n) from M{n) will be mapped to o* >(2 5{n) — {k — 1)) at the mirrored position in the second 
half of the later defined first fundamental unit {N = 1). To find this odd number the first difference 
set A A^(n), introduced in connection with Ai[n) from the definition 3, becomes important in order to 
count the gaps in the sequence of odd numbers when they are reduced modn. o* >(2 5(n) — {k — 1)) is 
found by adding to mk{n) twice the value of the sum of the gaps from the position k to the center, the 
mirror axis. This is 2 (Amfc_|_i(n) + Amfc+2(ri) + ... + Am5(„)(n) + 2). The 2 in this sum is half the 
gap-length from the value n — 2 at the position 5{n) and n + 2 at the next position, the mirror-position 
5{n) + 1 of 5{n). This is a telescopic sum which becomes 2(— mfc(n) + m^(„)(n) + 2). Therefore, the 
value of o* >(2(5(n) - (k - 1)) is + 2(-mfc(n) + m-5(„)(n) + 2) = 4 + 2m5(„)(n) - mk{n) = 
2n — mi^(n), because always nig^^-j{n) = n — 2 holds. See figure 5, n = 15 with the values 30 — 1 = 29, 
30 - 7 = 23, 30 - 11 = 19 and 30 - 13 = 17 for k = 1, 2, 3 and 4, respectively. Now it is 
clear that this mirroring leads to the correct number of reduced odd numbers for the second half of the 
fundamental = 1 unit, because gcd{mk{n),n) = 1, as member of 7W(n), implies for the mirror image 
also gcd{o^ >(2 5{n) — {k — 1)), n) = 1. All positive odd numbers relatively prime to n and not exceeding 
2n — mi{n) = 2n — 1 have thus be found. 

Proof of ii): From above we know that the number at position k = 2 5{n) + 1 is 2n + 1 because the 
one for k = 2 6{n) has been shown to be 2 n — 1 and gcd{2 n + 1, n) = 1 (indirect proof by assuming the 
contrary, using odd n; this is similar to the proof given in connection with definition 3 of A4{n)). Now it 
is clear that a shift in the o* > numbers with P(n) := 2n leads from the fundamental unit No. = 1 
to the second one, = 1, by putting o* >(2(5(n) + k) = o^^^{k) + P{n), for k G {1, 2, (5(n)}. 
This is obvious because the gcd value 1 is not changed by adding 2n. This process can be iterated 
to find the mod2n periodicity structure stated in ii). See figure 3, n = 15 with fc = 14, iV = 2: 
0*5^^(14) = 49 = 2 • 15 + >(14 - 1 • 2 • 4) = 30 + 0^5 > (6) = 30 + 19. □ 
The sequences {o* >} for n = 15 and n = 21 are found in A007775 and A206547 , respectively. 

The sequences of the Moddn residues of the numbers o* >(A;), for k = 1, 2, for prime moduli 
n = p = 3, 5, 7, 11, 13, 17, and for the first odd composed ones for n = 15, 21 are found in A000012 , 
A084101 , A110551 , A206543 . A206544 , A206545 , and A206546 . A206548 . respectively. 
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For the following formula for the nonnegative odd numbers reduced modn involving floor functions we 
need the following list (increasingly ordered set) J-{n) of length L{n) = 2 5{n) derived from the list AA^. 

Definition 5: List F{a) 

F{n) = {/(n,l),...,/(n,2 5(n)} with 

/(n,j) = forj G {l,2,...,<^(n)-l} 

f{n,6{n)) = 1, f{n,6{n) + j) = f{n,6{n) - j), for j G {1, 2, (5(n) - 1}, and 
fin,2 5{n)) = 0. (85) 

Here Amfc(n) := mk{n) — mfc_i(n), (the kih. element of AA^(n)). This list J-{n) is obtained from first 
enlarging AA4 by mirroring the first 5{n) entries at the last element 4, to obtain a list of order 2 5{n) + 1. 
Then the first and last element of this new list is put to zero and all the other elements are diminished 
by 2, thus obtaining a set of only even numbers. Then one divides by 2, omits the final 0, and reverses 
the remaining list. 

Example 16: -F(15) 

n = 15 = 3 • 5, 2 5(15) =2-4 = 8, A^(15) = {1, 7, 11, 13}, 7W(15) = {0, 1, 7, 11, 13, 17}, AA?(15) = 
{1, 6, 4, 2, 4}, the mirror extension is {1, 6, 4, 2, 4, 2, 4, 6, 1}, the reduction step leads to {0, 4, 2, 0, 2, 0, 2, 4, 0}, 
and finally, dividing by 2, omitting the last and reverting, leads to J^(15) = {2, 1, 0, 1,0, 1, 2, 0}. Except 
for the last there is a mirror-symmetry around the fourth entry 1. See the first row of table 5. 

For the odd squarefree composite numbers n =A024556(m), m = 1, 2, ...,17, see table 5 for A4{n), 
AM{n) and T{n). 

Proposition 13: O.g.f. and formula for O* > elements. 



oo 



With Gn{x) := V o;>(A;)a;^ for odd n, one has 



(l-3;2'5H)(l-x) 



2 5(n)-l 

1 + ^<>(^ + l)2;^' + J. (86) 

k=l 



k + jj- 1) 
2(^(n) 



(87) 



with the first differences Ao* >(j) := o* >(j) — o* >(j — 1) 
This generates 

2S{n) 

ol^{k) = 2fc-l + 2 fi^^j) 

i=i 

Note that the numerator polynomial in Gn{x) needs the first differences of the sequence members of the 
fundamental = 1 unit, which due to the mirror symmetry of proposition 12 can be reduced to the 
first differences of {7712(71), 7^3(71), .., 77i5(„),7i + 2} with 7)t,5(„) = n — 2. 

One could use, like in proposition 11, o„^>(0) = and let the sum in Gn{x) start with A; = 0. 

Proof: We start with the periodicity of the sequence due to proposition 12 ii). Taking the difference 
On,>{k) — o* >(fc — 1) = o'^^^{k — 25{n)) — o* >(/c — 1 — 2^(77)) leads to the recurrence (we omit 
all unnecessary indices) o{k) = o{k — 1) + o{k — 2 6{n)) — o{k — 1 — 2 5{n)), for k > 2 6{n) + 1. 
Here one needs also o* >(0) := —1. For some of these recurrences the o.g.f.s were determined by 

2 5(n) 00 

R. J. Mathar in e.g., A045572 and A162699 . I general Gn{x) = ^ o{k) x'' + x ^ o{k - 

k=l k=2S{n)+l 

00 00 

l)x'=-i + x^-'W Y 0(77 - 2(5(77)) - Y 0(77 - 2(5(77) - l)x"-2'5(")-i, where the 

fc=2 5(n)+l fc=2 5(n)+l 
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recurrence has been used (and the infinite sum has been reordered, not bothering about absolute con- 

2 5{n) 

vergence, in the sense of formal power series). Shifting the indices, one arrives at = o{k) + 

k=l 

X {g„(x) - (o(l)x + .... + o(2 5(n) - 1)x2^(")-i)} + x2^(")G„,(x)-x2'5(")+i (G„(x) + (-1)), where 



the —1 resulted from putting o(0) 



-1 (see above). This rearranges into (1 



,2S(n) 



+ 



x'^^(-^)+^)Gn{x) = X (1 + Ao(2)x + ... + Ao(2(5(n))x2'5(")-i + x^-^^")). Factorizing the bracket on 
the l.h.s. and division leads to the claimed form for Gn{x). Of course, the denominator factor (1 — x^'^'-'^-') 
can be factorized into cyclotomic polynomials. 

For the proof of the second part we use L := 2 6{n). Given G„(x) with the L + 1 input coefficients of the 
numerator polynomial one derives the claimed explicit form for o* >(A;), by defining first the sequence with 

~'~ , generated by 7 (partial sums of the characteristic sequence 



entries 6_l(A;) := 



L 



(1 



(1-x) 



for multiples of L, then shifted). The numerator polynomial leads for o* >(A;) to a sum of L + 1 floor- 



functions with decreasing arguments, starting with 
k-1 



k-O+L- 1 
Z 



ending with 



k-L+L-1 
I 



L 



and corresponding coefficients. In order to find a standard form for this sum we use the following 



Jioor-identity 



L-l 

E 

j=0 



k + j 
L 



L-l 



which can be seen from > {k + j) (modL) 



k, for L £ 



L(L-l) 



and k £ 



(88) 



j=0 



/c S Z, which is trivial (just add the L terms 



which are 0, 1, ... , L — 1 in a certain cyclic order), and the relation between modL and the floor function 
k 1 

k{modL)). This identity allows us to lower L consecutive coefficients of this sum of 



L 



L + 1 terms by 1, producing a term k if one uses the identity for the first L terms because the second to 
last floor-argument is then |_^^J , and the identity is read backwards. If one uses the identity for the last 
L terms one produces a fc — 1, because the last term has |_^^J . When we apply this identity twice in the 
described way we pick up n -|- (n — 1) = 2n — 1 and the first and last coefficient, which were originally 
1, become 0, and all other coefficients are diminished by 2 because they participate in both applications 
of the identity. It is guaranteed that all coefficients are now even and > 0, because the coefficients 
except the first and last one were even and > 2, because these numerator polynomial coefficients resulted 
from first differences of the sequence of odd numbers. Therefore one can extract a factor 2 and if the 

k 

floor-functions are written with increasing arguments, starting with 



L 



(the second to last term in the 



original order), one has exactly the coefficients given by the list J-{n) of length L of definition 5. This 
proves the explicit form for o* >(A;). □ 

Example 17: O.g.f. G8(x) 

{1 + 2{x + x'^ + x^ + x^ + x^ + x'^ + x^) + 



Gs{x) 



1 + X 



(89) 



(1 - x8) (1 - X) 

generating, with offset 0, the odd numbers A004273 . 

Before coming to the multiplicative group Moddn we define reduced residue systems Moddn as well as 
reduced odd residue systems Moddn. See e.g., [2], p. 113 for the modn case. 
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Definition 6a: Reduced residue system Modd n 

A reduced residue system Moddn {RRSn) is any set of ip{n) pairwise incongruent Moddn numbers, each 
of which is relatively prime to n. 

E.g., n = 15, 93(15) = 8, {1, 2, 4, 7, 8, 11, 13, 14} or {29, 32, 26, 23, 22, 41, 17, 16}, etc. The first one is the 
smallest positive one. These systems will not play a role later on. 

Definition 6b: Reduced odd residue system Modd n 

A reduced odd residue system Moddn (RoddRSn) is any set of 5{n) odd pairwise incongruent Moddn 
numbers, each of which is relatively prime to n 

E.g., n = 15, 5(15) = 4, {1,7,11,13} or {29,27,19,17}. The first one coincides with 7W(15) from 
eq. (65) and is the smallest positive reduced residue system Moddn. Remember that the parity of the 
members in each Moddn residue class, not in class [0], is the same (see lemma 12). Later on we will 
restrict ourselves mostly to the system Ai (n) . 

We next study the multiplicative group Moddn. The elements are the residue classes [nrij] , for j = 
1, 2, 5{n), corresponding to the members of the reduced residue system A^(n). In fact, we will take 
these representatives, multiplying Moddn, as we have done above. 

In order to see the group structure one first convinces oneself that this set M{n) is closed under Moddn 
multiplication. This follows from proposition 10 and lemma 12 which showed that these classes have 
only odd numbers, and every odd numbers appear exactly once because of the definition of these classes. 
The associativity, commutativity and the identity element 1 are also clear. We do not have a formula 
how to find the inverse element ni~^ of nij but coming to this Moddn multiplication from the study of 
automorphisms of the splitting field for the minimal polynomial C{n,x) for the algebraic number p{n), 
it is clear that these inverses have to exist (the invariance property defines a group). We are dealing here 
with the Galois groups for these polynomials. The cycle structure for n = 1, 2, 40 is given in table 
6. Because each group is Abelian one has for every group of prime order a cyclic group, which can be 
checked for the examples given in table 6. 

Remark 6: The multiplicative group Moddn is cyclic if 5{n) is prime. 

This is a corollary on the fundamental theorem on finite Abelian groups (see e.g., [28], p. 49, see also 
[6], p. 511, Cauchys Theorem A. 1.5), or use Lagrange's theorem on the order of subgroups, and the fact 
that the powers of an element, not the identity, generate a cyclic subgroup, to prove that in fact every 
group of prime order is cyclic. The values n for which the order 6{n) is prime are given in A215046 . 
This is the sequence [4,5,6,7,9,11,23,47,59,83,107,167,179,227,263,347,359,...]. Of course, there are 
other values n with cyclic Moddn group, like n = 2, 3, 8, 10, .... For n = 1 one has the trivial case with 
cycle structure [[0]], also a cyclic group, viz Z\. One can give a more general sequence of n numbers 
with cyclic Qn — multiplicative Moddn group, namely A210845 . 

Remark 7: The multiplicative group Moddn is cyclic if 6{n) is squarefree. 

The squarefree numbers are given in A005117 , namely [1, 2, 3, 4, 5, 6, 7, 9, 11, 13, 14, 18, 21, 22, 23, 25, 29, 
31, ...]. This remark 7 results from the fact that A000688 , giving the number of Abelian groups of order 
n, is 1 exactly for the squarefree numbers A005117 . See the formula, based on the H.-E. Richert reference 
quoted there. Because for each order there is at least the cyclic group these values lead necessarily to a 
cyclic group. The above given A215046 values are a proper subset of those from A210845 . There are. 
however, still more values n with a cyclic Moddn group. Missing are e.g., 8, 10, 15, 16, 17, .... All the n 
values with cyclic Modd n group are in A206551 . The complementary sequence is A206552 , giving the n 
values with non-cyclic Moddn group. See Table 8 for all values n < 100. 

As an aside we remark that special squarefree numbers are the so called cyclic numbers A003277 . It is 
known that if the order of a (finite) group is a cyclic number then there is only one group, the cyclic one. 
See a comment on A003277 . Therefore AOOOOOl (n), the number of groups of order n is 1 if n is a cyclic 
number, in fact, the reverse also holds: if AOOOOOl(n) = 1 then n is a cyclic number. This can be taken 
as an alternative definition for cyclic numbers because then there is only one (non-isomorphic) group of 



34 



this order which has to be the cychc group. See also Yimin Ge's Math Blog [YiminGe], where the 'only 
if statement in the proposition may be misleading but in the proof the given statement is correct. E.g., 
for order 6 (not a cyclic number) there are groups other than Zq. In fact A000001(6) = 2, and there 
is the (non-Abelian) group D(3) (dihedral group). (It is clear that cyclic numbers are not the numbers 
n for which the multiplicative group modn, which is Abelian of order <p{n), is cyclic. These numbers 
are given in A033948). Because we are dealing with Abelian groups the squarefree numbers are more 
interesting here. 

Recall the situation for the Galois group for the cyclotomic polynomials (the minimal polynomials for 
Cin), an n-th root of unity) which is isomorphic to the multiplicative (Abelian) group modulo n. See Table 
7 were we have listed these non-cyclic groups for n G {1, ... , 100}. This n— value sequence is known as 
A033949 . (The case n = 15 is an exercise in [7], p. 159, Ex. 9.6. 3)a)). The values n with a cyclic Galois 
group coincide with the moduli n which possess primitive roots r = r{n), i.e., the order of r modulo n 
is if. r"^^^^ = l[modn), and no smaller positive exponent k satisfies this congruence. These moduli n 
are known to be exactly p'^, 2pf , 1 and 2, with some positive powers e and e'. See e.g., [23], Theorem 
2.4-1, p. 104. All other n lead to non-cyclic multiplicative groups modulo n . See A033949 . For the 
smallest primitive roots in this case see [1], pp. 864-869, the column called g. We do not know a similar 
characterization of the non-cyclic numbers n for the M)ddn group which are shown in Table 8. However, 
these numbers have also to appear in A033949 . It is clear that we should determine primitive roots 
r = r{n) for the Mjddn multiplication, i.e., find those r from Ai{n) which have order 5{n): the smallest 
positive /c which satisfies r*^ = +l{M)ddn)isk = 6{n). See A206550 for these smallest positive primitive 
roots Moddn, where a entry, except for n = 1, indicates that there exists no such primitive root. This 
sequence starts, with offset 1, as [0, 1, 1, 3, 3, 5, 3, 3, 5, 3, 3, 0, 7, 5, 7, 3, 3, 5, 3, 0, 11, 3, 3, 0, 3, 7, ...]. In general 
one does not expect a formula for these primitive roots r = r{n), because also in the multiplication 
modulo n case there is no one available. 

Proposition 14: Number of Modd n primitive roots 

If a primitive root r for the multiplicative group Moddn exists there are ip{5{n)) of them. 

This is the sequence A216322 . For example, n = 13, 5(13) = 6, 97(6) = 2 with ri = 7 and r2 = 11; 
n = 14, (5(14) = 6, ^{Q) = 2 with n = 5 and ra = 11. 

Proof: If a M)ddn primitive root exists then the multiplicative group Mjddn is cyclic. The order of this 
group is d{n), Then the number of pairwise incongruent primitive roots is obtained exactly like in the 
case of primitive m-roots of unity on the unit circle by the number of relatively prime numbers less than 
m which is Eulefs ip{m). Here m = 6{n). □ 

Now the question whether the multiplicative group Moddp with p a prime is cyclic is answered. Up to 

p - 1 

now we know from remark 7 the positive answer only for those primes with 6{p) = squarefree. 

These are the primes given in A066651 . The case p = 2 with 5 (2) = 1 is trivially cyclic. For general 
prime p one wants to show that there exists a primitive root in the multiplicative Moddp group. This is 
analog to the modp case, where a proof can be found in [23], Theorem 2.36, p. 99, or in [2], ch. 10.4, 
pp 206 ff. One has, however, to be careful to use only the multiplicative group structure in the proof. 
Already in connection with lemma 10 we have given a warning to use theorems valid in the modular 
arithmetic modn in our case. Another failure occurs for the theorem [23], Theorem 2.3(c), p. 49, e.g., 

X = y {mod nil) , X = y{modm2) if and only if x = ?/ (mod — — - — - ) . For the Moddn case 



gcd{mi,m2) 

one has the counterexample 19 = 5 (Mode? 3) and 19 = 5 (Modd 7)) but 19 is obviously not congruent 
5 (Modd 21). Several lemmata, the analoga of mod n facts, are collected before stating proposition 15. 

Lemma 19: Analogon of the Fermat-Euler Theorem 

If a is odd and gcd{a,n) = 1 then a*^^"-* = +1 (Moddn), with delta{n) = A055034(n) (see the start of 
section 3). 
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Proof: Analog to e.g., [2], Theorem 5.17, p. 113, with m ^ n, <p ^ 6, the order of the set (reduced 
odd residue system) Ai{n) of eq. (65), and modulo — t- Moddn (we avoid the term Moddulo). Instead 
of M{n) one can use any other reduced odd residue system Moddn (see definition 6b), with mj{n), j G 
{1, 2, ... , S{n)}, replaced by any member of the residue class [mj{n)]. The cancellation used at the end of 
the proof in [2] can be done here by multiplying with the existing inverses of the biS there. We are dealing 
with the multiplicative group Modd n. 

Next follows the definition of the order of an element from a reduced odd residue system Moddn. 
Definition 7: Modd n order of a from a RoddRSn 

The Modd n order of a positive (odd) integer a from a reduced odd residue system Modd n is the smallest 
positive integer h such that = 1 (Moddn). 

E.g., n = W,a = 9,h = 2 because 9^ = 9 (Modd 10), and 9^ = 1 (Modd 10). See the table A216320 
corresponding to the Moddn orders h of the smallest positive RoddRSn members A216319 . 

Lemma 19 guarantees the existence of a Moddn order h < 6{n) for each positive odd number a with 
gcd{a,n) = 1. 

Example 18: n = 12, 6(12) = 4, A^(12) = {1,5,7,11}, 1^ = 1 (Modd 12), 5^ = 1 (Modd 12) ^ 5^ = 
(Modd 12), similarly for 7 and 11. See the cycle structure for n = 12 in table 6. 

Lemma 20: The Moddn order divides 8{n) 

If h is the Moddn order of a then h\8{n). Moreover, = [Moddn), w.l.o.g. j > k, if and only if 
h\{j-k). 

Proof: The first part is analog to [23], Corollary 2.32, p. 98, with ip ^ 5. It uses lemma 19, the 
Fermat-Euler analogon. Remember that always gcd{a, n) = 1 from the order definition. The second 
part, for which one can use the group property, follows from the analog of [23], Lemma 2.31, p. 98. (In 
the older German version of this book, vol. I, the analog of this lemma 20 appears as Satz 22.3 on p. 63.) 
□ 

Lemma 21: On the Moddp order of a, p a prime 

If h is the Moddp order of a, with p a prime, one has (a'^)'* = 1 (Moddp) for all k, and 1 = 
a^, , a?, ... ,a^~^ are pairwise incongruent Moddp. 

Proof: First part: (a'^)^ = (a^)^, = 1 (Moddp) and Moddp respects multiplication (see proposition 
10). Second part: assume the contrary, i.e., a* = (Moddp), l<j<i<h — 1. Apply lemma 20 ior 
n = p, showing that /i|(i—j), but i — j <i< /i — 1 leading to a contradiction. □ 

Lemma 22: X*^ = 1 (Moddp) has at most h incongruent solutions 

The congruence = 1 (Moddp), h a positive integer, has at most h pairwise incongruent solutions 
Moddp. 

X^ 

Proof: If is even then X^ = l(modp). and due to the modp theorem, e.g., [21], Theorem 



p 

42, p. 80, there are at most h incongruent solutions. If 



is odd then —X = l(modp), i.e.. 



P 

X^ = (p—1) (modp) which again has at most /i incongruent solutions. □ 

This is a weak statement, but sufficient for the following. One could try to prove that the number of 

p — 1 

Modd n incongruent solutions of the congruence X^ = 1 (Moddp) is gcd(h, — - — ) if p is odd. This 

number is trivially 1 for p = 2. 

Lemma 23: On the Modd n order of a*^ 

If h is the Moddn order of a, then — — - — — is the Modd n order of . 

gcd(h, k) 

Proof: This is the analogon of [23] lemma 2.33, p. 98, using the present lemma 20, part 2. 
Now we are ready for the following proposition. 
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Proposition 15: Existence of Modd p primitive roots 

For every prime p there exists a primitive root for the multiplicative group Modd p. 

Proof: One shows that there are precisely (p{5{p)) primitive roots and because this number is always 
> 1 the claim will then follow. 

Following the proof of [23], Theorem 2.36, p. 99, we infer from the present lemma 19 (with n ^ p) that 

p-l 

every positive odd integer a with gcd{a,p) = 1 has a Moddp order h = hp{a), 1 < h < d{p) = — - — , 

and from lemma 20 h\5{p). Moreover, {a^)^ = \{Modd p) for all k, and 1, a^, ...a^~^ are h pairwise 
incongruent odd numbers Moddp. From lemma 22 these are all the solutions of the congruence = 1 
{Moddp). Lemma 23 shows that there are '^{h) numbers , with k G {1, ... /i — 1} which have Moddp 
order h = hp{a) because f{h) is the number of such k with gcd{k,h) = 1. The case /i = 1 is covered 
with f{l) := 1. 



With [2], p. 207, we define for each h\5{p) the set Ap{h) 



a S J^{p) and Moddp order of a is /i 



There are t{5{p)) (r = A000005 , the number of divisors) such disjoint sets. This is the sequence 
[1, 1, 2, 2, 2, 4, 4, 3, 2, 4, 4, 6, ...] = A216326. E.g., p = 7, M{7) = {1, 3, 5}, r(3) = 2 (see also 
row No. n =p = 7 of the Moddn order table A21630). Call ipp{h) := |Ap(/i)|, the number of elements 
of this set. We have seen that ipp{h) is or ip{h) because for each h\5{p) the set Ap{h) is either empty 
(no a € M{p) has Moddp order h) or this set has <^(/i) elements. Thus ipp{h) < (p{h). It is clear that 
ipp{h) = 5{p) because each a G M.{p) belongs to one of these sets Ap{h). Now a standard result is 

h\5(h) 

^p{h) = 5{p) {see e.g., [2], Theorem 2.2, p. 26, wit hn — )• 5(p)), therefore, {ip{h) — ■ipp{h)) = 0, 

h\5{h) h\S{h) 

but because (p{h) — ipp{h) > it follows that ^p{h) = (p{h) > 1, not 0, for each h\6{p). This holds 
especially for h = 6{p). Hence the number of primitive M)ddp roots is ipp{5{p)) = ip{5{p)) > 1, and the 
proof is complete. □ 

For the sequence 5{p{n))n > 1, see A130290 . For the smallest positive primitive Moddn roots see the 
sequence A206550 . Here for prime n. 

Corollary 8: Cyclic Modd p group 

The multiplicative (Abelian) Moddp group, for p a prime, is the cyclic group .^^(p). 

We now concentrate on those groups Moddn which are cyclic and notice that whenever 6{n) is even there 
exists a unique smallest positive odd number > 1 which solves the congruence = 1 {Moddn). This is 
the nontrivial solution of this congruence. The trivial one is x = 1 (standing for the class of solutions 
1 {Moddn) which also includes —1). We prove this first for even n. 

Proposition 16: Nontrivial square-root of 1 (Modd 2 k), k > 2 

If n = 2k, for k G {2,3,4,...}, and if the Galois group Gal{Q{p{n))/Q) is cyclic, a unique smallest 
positive solution r > 1 of the congruence = 1 {Modd2 k) exists, and it is r = n — 1. 

Proof: For even n = 2 k, k > 2, the cyclic group is then Zs(^2k) and 6{2k) is even because with the 

standard prime factorization (5(2^ • • -p^j^) = 2^^"^ p^ ^ {pj — 1), and if = then e > 2, and 

j=l..N 

if > 1 then e > 1. In both cases d is even, because in the second case there is at least one j and 
Pj — 1 is even. Even order 6{2 k) of the cyclic group means that powers of the generator c of the J-cycle, 
the smallest positive primitive root (of 1) (Modd 2 k), come in pairs which are mutually inverse, except 

for 1 = c'^(^'^) and s := c 2 . This s > 1 is unique and solves of the congruence = 1 {Modd2 k). 
(n - 1)2 

Because = n — 2 + is, with even n, even, and in order to compute (n — 1)^ {Moddn) 



n 
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one has to determine (n — 1)^ (modn) , which is 1. Therefore, the unique s for even n is s = n — 1. 
□ 

Example 19: Nontrivial square-root of l(Modd8) 

49 



6, hence P8(49) = +1, and 49 = +l{mod8), 



s{8) = 7: 72 = 49 = l{Modd8), because 
therefore 49 = +l{Modd8). 

p — 1 p — 1 
Similarly, for n = p, p an odd prime, the cycle length is 6{p) = . In the event that is even, 

i.e., Hp is of the form 4/c + l, i.e., p = 1 (mod 4), (see A002144) a smallest positive nontrivial solution 
of the quadratic congruence = 1 (Moddn) exists. We will numerate the primes congruent 1 {mod 4) 
by defining p{n) := A002144(n), n > 1. E.g., for n = p{2) = 13 one has, besides the trivial solution 
X = 1 (ModdlS) (this residue class includes —1), x = 5 (ModdlS) as a solution. This means that in the 
group table the row for the element 5 has a 1 in the diagonal, like in the row for the identity element 1. 
These elements are self-inverse like the identity element 1. Note that in the ordinary reduced modn case 
there is always a 1 in the diagonal for row p — 1, besides the 1 in the row for the identity element, but in 
that case p—1 = —1 (modn), and hence is a trivial solution. See A206549 for these non-trivial solutions 
in our case. The other inverses in these cyclic groups can be read off the cycle structure. One pairs the 
first entry, the generator c of the cycle, with the second to last (the one before the 1, the second with 
the third to last, etc. E.g., for n = 13 one has the inverses 7~^ = 11 (Modd 13), 11"^ = 7 (Mode? 13), 
and 3~^ = 9 {M)dd 13) , 9~^ = 3 (Modd 13), and the left over 5 is self-inverse, as mentioned above. This 
pairing of inverses is clear from the cyclic structure. 

The next objective is to find some method to compute this uniquely existing smallest positive non- 
trivial solution s{p) of the congruence x^ = \ (Moddp) for given prime p = p{n) = A002144(n) for 
the k{n) := Eil}}^ values given in A005098 . This A;-sequence starts with [1,3,4,7,9, 10, 13, 15, 18, ...]. 
The following Z-algorithm will determine all these non-trivial solutions s{p{n)) of the congruence x^ = 
1 [M oddp{n)). We start with the Ansatz s{p) = y^(o + l)p — 1, with some odd number o. This 

+ 1 

means that we are looking for primes of the form , and we are interested only in primes of the 

+ 1 

form Ak + 1, i.e., p{n). Because s(p)^ = (o + l)p — 1, s{p) is odd, say 2K + 1. This implies 

2K (K + 1) + 1 = — — p, which shows that has to be odd, i.e., o = 4/ + 1. Now we compute 

v — 1 (2 K + l)"^ — (Al + Dv 
the even number 2k = '-^ = ' ^ ^ '— = 4T{K) - 21 - 2k{Al + 1), where 

we have introduced the triangular numbers T{K) := ■^(^+^) ^ given in A000217 . This implies that 
2k{2l + 1) = 2T{K) — I . Therefore, I is even, say 21, and we have, with inputs k, I and K: 
p = Ak + 1, s^(j5) = (0+ l)p — I, o = 81 + 1 and s{p) = 2K{p) + 1 . Because now 
2T{K) = 2k + 2l{Ak + 1) = 21 + 2k{Al + 1), we find the 'nice equation' 

AT{K) + I = {Ak + 1){AI + 1) . (90) 

For given prime p of the form Ak + 1, i.e. for k = ^ (see A005098), we determine the minimal / G No, 
such that A; (1 + 4/) + / produces a triangular number, namely T{K) . The index of this triangular 

number is then K = V^'^^^^^~ tJ_ — 1^ have uses Maple [18] to compute these minimal I values. It 
will be seen from the later proposition 15 that there is always such a minimal /, and not all non-negative 
values appear. Precisely the entries of A094178 and occur as minimal I values. We will from now on 
denote this minimal I also by I (hoping to cause no confusion). If Z = then A; is a triangular number, 
and conversely. The non-vanishing / values are characterized by having as elements of the squarefree 
kernel set of 4 / + 1 only distinct primes congruent 1 (mod 4). E.g., 1 = 2 will not appear because the 
squarefree kernel of 9 consists of the prime 3 which is not congruent 1 (mod 4) . In order to proof that 
only numbers from A094178 can appear as positive I values, we use the 'nice equation' and observe that 
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4:T(K) + 1 = + (K + 1)^, i.e., it is the sum to two neighboring squares with gcd{K,K + 1) = 1 
(shown by indirect proof, see the remark following definition 3 ). Now a theorem ensures that all prime 
factors of 4T(K) + 1 are congruent l(mo(i4). See the theorem 3.20 on p. 164 of the reference [23]. 
From the 'nice equation' we have either 41 + 1 = 1, i.e., I = 0, or, if I is non-vanishing, then the 
squarefree kernel set of 4Z + 1 has only primes congruent 1 (mod 4). Therefore I can be or it is from 
A094178 , i.e., from [0, 1, 3, 4, 6, 7, 9, 10, 13, 15, 16, 18, 21, 22, 24, 25, 27, 28, ...]. The / = value appears 
for the primes [5, 13,41,61, 113, 181,313,421, ...], which is the sequence A027862 . In this case k = 2zl 
is a triangular number. / = 1 for the primes [17, 29, 53, 73, 109, 137, 281, 397, 449, 593, 757, ...] This is the 
sequence A207337(?i), n > 2. These are the primes p such that 5p = 4:T{K) + 1 for some K = K(p). 
E.g., 5 17 = 85 = 421 = l,thusir(17) = 6. These are also the primes of the form (m^+l)/10. Just use 

m = m{K) = 2 K + 1. In general, k = ^ — —, (see eq. (90) and K = ^ — — . 

E.g., p = 53 belongs to / = 1 and K = 11, and T[K) = 66. This checks with the 'nice equation'. We 
do not give here the tables of the sequences l{n) = l{p{n)) and K{n) = K{l{n)) (this incorrect renaming 
of arguments should not lead to confusion) for given p[n) because after the next proposition 17 we will 
find another way to obtain the K(n) numbers directly, hence the corresponding l(n) numbers, and from 
these the desired s{p{n)) solutions from s{p{n)) = 2K{n) + 1. 
Proposition 17: Congruence 4T(X) + 1 = + (X + 1)^ = O(modp) 

There are precisely two incongruent solutions of the congruence f{X) := 2X'^ + 2X + I = 0{modp), 
provided p = p{n) = A002144(n), n > 1. The smallest positive representative will be called K{n), and 
the next larger incongruent one is then K2{n) := p{n) — 1 — K[n). 

Proof: Because of the degree 2 of / this congruence has at most two incongruent solutions. We shall 
see that in fact there are two. This congruence is reduced to a problem of quadratic residues, following 
a standard prescription (see, e.g., [21], pp. 132-3). The discriminant of / is D = 2^ — 421 = —4. 
Multiplying / by 8 yields {AX + 2f + A = (modSp). With Y = AX + 2 this is Y"^ = D (modn), 
with a composite modulus n = 2^ p. This is a quadratic residue problem, but gcd{D,n) is not 1, but 4. 
Theorem 77 of [21] is applied with d = A, e = 2, f = 1, a+1 = —1 and ni = 2p Thus the problem is 
reduced to the quadratic residue problem = —1 {mod 2 p) . This is solved by studying the congruences 
for the powers of primes, here just 2 and p, separately (see e.g., [21], sect. 26, pp. 83-5). The congruence 
modulo 2 has only the solution +1 (because —1 = +l{mod2,)), and modulo p one can consult the 

Legendre symbol ^ — ^ which is (— 1)^2~ (see e.g., [23], Theorem 3.2(1), p. 132). Therefore —1 is a 

quadratic residue modulo p if and only if this symbol is +1, demanding that p ^ 1 {mod A), i.e., p = p 
from A002144 . Call the smallest positive solution xq, then p — xq is also an incongruent solution modulo 
p, and two is the maximal number of solutions because of the degree 2 of this congruence. This implies 
that there are 1-2 = 2 incongruent solutions of this congruence modulo 2p (see e.g., [21], Theorem 46, 
p. 84). Returning to the original problem this proves that there are also two incongruent solutions. If the 
smallest positive solution for p{n) = A207337(n) is called K{n), then the next larger incongruent solution 
oiAT{X) + I = X"^ + {X + 1)2) = 0{modp{n)) is K2{n) := p{n) - 1 - K{n), which is obvious. 
□ 

The pair of sequences of all positive solutions modulo 5, 13, and 17 are given in A047219 , A212160 and 
A212161 , respectively, where in each case the even indexed members are the positive solutions congruent 
to K{n) and the odd indexed ones are the positive solutions congruent to K2{n), i.e., a{2 k) = kp{n) + 
K{n) and a{2k + 1) = kp{n) + K2{n), /c > 0. For the three given examples n = 1, 2, 3 respectively. 
The companions {K{n), K2{n)) have been computed with Maple [18] and they can be found as A212353 
and A212354. The first entries are for /sT: [1,2,6,8,15,4,11,5,13,27,37,45,16,7,18,52,64,...] and for 
K2: [3, 10, 10, 20, 21, 36, 41, 55, 59, 61, 59, 55, 92, 105, 118, 96, 92, 126, ...] corresponding to the primes p: 
[5, 13, 17, 29, 37, 41, 53, 61, 73, 89, 97, 101, 109, 113, 137, 149, 157, 173, ...]. Note that proposition 17 yields 
directly the searched for nontrivial solution of s{p{n))'^ = +1 {Modd p{n)) via s{p{n)) = 2K{n) + 1. 
Compare this with the s values given in A206549 with first entries [3, 5, 13, 17, 31, 9, 23, 11, 27, 55, 75, 91, 33, 15, 37, IC 
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Because s(n) := s{p{n)) = 2K{n)+l < p{n)-2,i.e., K{n) < ^ ^ ^ iff 52(n) := K2{n) + l > p{n) 

and s(n) < s2(n), n > 1. Tlius only s(n) G tlie restricted odd residue class Modd p, and the 

solution s2{n) is discarded. 

This proposition and the 'nice equation' eq. (90) show that the / — algorithm from above will indeed pro- 
duce a solution /, related to the existing K{n) for given j5(n), via l{n) = ( ^ — ^ — — 1 ) /4. The 

V p(n) ; 

smallest positive non-trivial solution s(p(n)) of the congruence = -|-1 {Moddp{n)) is then 2K{n) + 1. 



4120, + 1 
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l)/4 



s{n) = 31, with 



25 = o and 31^ = 961 = -1 (mod 37) = 36. 



See K{n) = A212353(n), n > 1. E.g., n = 5, p{5) = 37, K{5) = 15, / 

312 
37 

As a application involving s{p{n)) we derive the analog of Wj'ison's theorem Yl T^ip) = {p ~ 1)! = 
{p— 1) (modp) = —1 (modp), for each prime p (see e.g., [2], Theorem 5.24, p. 116, or [10], Theorem 80, 
p. 68), where TZ{n)is the set of the representatives of the smallest positive reduced residue system modn 
for n > 2 which is of order (p{n). In the Moddp case we have to replace the set TZ{n) by the set M{n) of 
eq. (65) with order 6{n). 

Proposition 18: Analog of Wilson's theorem for Modd p 

{1 if is odd 
{Moddp), (91) 
s{n) if ig even, p = p{n), n G N 

where s{n) = s{p{n)), p{n) = A002144(n), stands for the above treated nontrivial root Moddn, i.e., the 
solution of the congruence = 1 [Moddn) which is not 1, if it exists. Note that —1 = -|-1 (Moddn), 
n > 2 (see lemma 10). For the double factorials [p — 2)!! see A207332 . 

Proof: The multiplicative group group Moddp, p a prime, is the cyclic group from proposition 

14. The order 5{p) of this group is even if and only Hp = p, i.e., a prime 1 (mod4). In this case we 
know from above that the non-trivial s{n) = s{p{n)) exists, and an algorithm for finding it has been 
given. Besides the unit element 1 and this s(n) all other factors in Y\ MiP) can be paired such that 
their product is 1 (Moddp) (see the discussion after proposition 16), leaving only s(p). In the other case, 
when p = 3 (modp), the group order is odd. Then all numbers besides 1 can be paired in the product to 
produce 1 (Moddp), and the result is therefore 1. □ 

For the cyclic Galois groups belonging to C(n,x) the list of the smallest positive primitive roots r(n) 
(i.e., the smallest element from M(n) which generate these cyclic groups) are found under A206550 . As 
mentioned above, we do not have a formula for those n > 3 values whose cycle does not start with 3 as 
smallest positive primitive root, like for n = [6, 9, 13, 14, 15, 18, 21, 26, 27, 33, 37, 38, 39, ...]. 

Open problems: 

• Proof of the conjecture on the g— sequence related to the discriminant of the minimal C(n,x) polyno- 
mials; or find a counterexample. 

• Characterization of the values n for which the Galois group Qn is non-cyclic. 

• Characterization of the values n for which the cycle of the cyclic multiplicative group Modd n = ^„ is 
generated by 3. 

• More theorems on multiplicative Moddn arithmetic. 
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Appendix A 



The proof of eq. (30), with replaced by the product-formula given there a bit later, is a (nice) 
application of PIE (the principle of inclusion and exclusion. See e.g., [5], Theorem 4.2, pp. 134 ff). 
As mentioned above, this formula uses only the distinct prime factors of n, the elements of the set 
sqfkset{n), the set of primes of the squarefree kernel of n. The product formula for can be read 
as the generating function for the elementary symmetric functions (here polynomials) in the variables 
= 1,2, ...,M{n). M{n), also called uj{n), is given in [27] as A001221(n). 



M{n) , , M(n) 



n 



n (i- -) = E(-ir-M-'-'^ 



(92) 



with (To = 1, and symbolically Ur = , where the sum extends over the ( ) terms with r 

^ p ... p J 

factors ^ with increasing indices. E.g., M = 3 with a2 = ^ + ^ + 

n-1 

To calculate k one starts, at the zeroth step (r = 0), with the unrestricted sum which is ^ (n — 1), 



fc=i 

gcd(fc,n) — 1 



and subtracts, in the next step (r =1), the sum of all multiples of each prime dividing n which are 

JL _ 1 

/ \ 

k = — I — — 1 ] , for j = 1, 2, M{n). This 
k=i ^ ^ 

leads, in step r = 1 to the subtraction of — ( 1 | . Now in this subtraction all multiples of the 

2 j \Pj J 

product of two different pjS appeared twice, therefore one has, in the next step (r = 2), to add them 

— 1 

Pi f J / X 

once. This is done by +1 ViVj ^ ~ '2 ( ~ step r = 3 one concentrates 

i<j k=l i<j ^P P- ^ 

on products of three different primes PiiPi2Pi3 with ii < i2 < ^3. Now such a 3— product appeared 
once in step r = ( in the unrestricted sum), —3 times in step r = 1, originating from the multiples 
{Pii Pi2 )-Pia^ (Pii Pis ) ■ Pi2 and {pi^ Pis) 'Pin and +3 times in step r — 2, from the multiples pi^ ■ {pi2 Pi-^), 
Pii ■ {PiiPiz)-: and Pig • {pij^Pi^)- Therefore, up to this stage, each such 3— product appeared once too 
much, and it is subtracted in this step r = 3 when all these 3— products terms are summed. Now the 
pattern starts to become clear. Up to, and including, step r = 3 one has for each 4— product the counting 
1— 4 + 6 — 4 = —1. Hence in step r = 4 one adds once the sum over all multiples of each such 
4— product. In general, in step r this counting will produce (—1)''+^ (from the alternating row r = 4 in 
the Pascal triangle A007318), and one will therefore add (—1)*^ times the sum over all multiples of each 
such r— product. This yields 



n-1 

gcd{k,n) — l 



Pii ) .^^.^ \Pi1Pi2 



il<...<iM(n) 



n 



P1P2 ■■■ PM{n) 



(93) 



Now the ±1 terms all cancel if one takes into account the number of terms of each r— sum, which is (f). 
This is due to the vanishing alternating sum over Pascal triangle's row M. Thus the result becomes the 
alternating sum over the elementary symmetric polynomials with the reciprocal distinct prime factors of 
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n as variables. 

n-l _2 ^(") 

(94) 



n-i 2 ^"^"^ / I 1 \ 



gcd{fc ,ri) — 1 

As mentioned above this is written as the generating function of these polynomials and the final result 
is with the second part of eq. (92) 

n-l 2 / -I \ 

s * - T n - ■ <'^> 

gc(i(A;,n) = l 

This shows that s(?i) := — /c is indeed , which appears in [27] as A076512 (n) / A109395 (ra). 

^ t,— 1 ™ 

pcd(fc,7i) — 1 

s(n) is identical for all n with the same distinct prime factors, independent of their multiplicity, i.e., it 
depends only on the elements of sqfset{n). Prom the multiplicativity of if one sees that this scaled sum 
s{n) is also multiplicative. E.g., s{12) = s(6) = s{2) s(3). 

Appendix B 

For the proof of eq. (56) we follow [24], p. 33, and consider the rewritten version 

71-1 ( 1 if n is odd , 

n (l - e^-'^) = (96) 
k=i I if 77, is even . 



In order to see that this is identical to eq. (56), just extract e'^* leading to the (2 cos) factors under 
the product, and in front summing in the exponent leads to the factor gz'^*^""^^ = z""^ = (—1) 
The l.h.s. of eq. (96) is then computed using (for the last step see e.g., [9], Exercise 50, p. 149) 



2 



1 + z + + ... + z 



n-l _ ^ ^ 



^ = n(-- e'-'-), (97) 



k=l 

n-l 



where one specializes to z = — 1. (By putting z = +1 one finds the result 2 sin I — | 

k=i \ 



n , n >2, 



with 1 for n = 1.) 
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Table 1: Reduced DSR-algebras (over Q), n = 3, 12. 



n 



reduced DSR-algebra 



S{n) 



DSR - basis 



10 



11 



12 



^(1 + ^) 



V2 



2 cos (f J 



2cos(f) 



(py/3 - (p 



2 cos (^) 



^2 + \/3 



-- 

pa 
{a - 

P' - 
(a - 

P' - 
pa 

(a - 

[CO - 

P' - 
r2 = 

pa 

a T 

TiO 

{a - 

P' - 
pa 

[u - 

{a - 



1 

2 

P + 1 

3, [a 



2] 



1+cr, a^ = I + p + a, pa = p + a 

1 + cj, cr2 = 2(T + 1, = 2(1 + a), 
■ p + T, pT = 2 a, a T = 2 p + T, 
1 + V2, T = V2p) 



I + a, a^ = 2 + p + a, pa = 2p + 1, 
p2 - 1), [r = p{a - 1) = 1 + p] 



I + a, a^ 



-1 + 3a, 



-1 + 4(7, 



-- p + T, pT = 3a — 2, ar = p + 2t, 
p^ - I = 1 + ^ , T = p{a - 1), 
(1 + p) V3^) 
2(-l + a) = 2^] 

1+cr, fJ^ = 1 + cr + W, 

l+u+r+w, uj'^ = l+p+a+T+uj, 
- p + T, pr = a + oj, puj = T + uj, 
-- p + T + oj, a oj = a + T + UJ, 
-- p + a + T + ijj, 

p^ - I, T = p{a - I), UJ = a{a - I) - I) 

I + a, ^2 = 2(1 + r), t2 = 3(1 + a), 
-- p + T, pT = 1 + 2a, ar = 3 p + T, 
1 + a, X = 2p], 

V2 



1 = 1 + V3, T 



(3 + ^3)) 



< 1 > 

< l,/0 > 

< l,/0 > 

< l,/0 > 

<l,p,a > 
<l,p,a,T> 

<l,p,a> 
<l,p,a,T > 



< l,p,a,T,uj > 



<l,p,a,T > 



n 



p(n) := 2 cos(-) , = S{k - I, p{n)), k=l,.... 



Ri = 1, i?2 = P, Rs = a, Ri = T, R^ = UJ, Rq = X (dependence on n suppressed). 

In round brackets the values for the basis elements are given in terms of p. In square brackets the linear 
dependent DSRs are given. Boxed n-numbers indicate linear dependent DSRs. 
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Table 2: Minimal polynomials of 2 cos ( — ) for n = 1,2,. ..,30. 



n 


C(n,x) 


1 


x + 2 


2 


X 


3 


x-l 


4 


x^ -2 


5 


X^ — X — 1 


6 


x^-3 


7 


x"^ — x'^ — 2x + 1 


8 


x^-4x^ + 2 


9 


- 3x - 1 


10 


x^ — 5 x^ + 5 


11 


x^ - x^ - 4 x^ + 3 x^ + 3 X - 1 


12 


x^ - 4x^ + 1 


13 


x*^ - x^ - 5 x^ + 4 x^ + 6 x^ - 3 X - 1 


14 


x^ - 7x^ + 14x^ - 7 


15 


x"' + x^ — 4 x^ — 4 X + 1 


16 


x« - Sx*^ + 20x4 - 16x^ + 2 


17 


x^ - x'^ - 7 x^ + 6 x^ + 15 x'' - 10 x^ - 10 x^ + 4 X + 1 


18 


x^ - ex"^ + 9x^ - 3 


19 


x^ - X** - Sx''' + 7x^ + 21x^ - ISx"^ - 20x^ + lOx^ + 5x - 1 


20 


x^ - Sx*^ + igx'* - 12x^ + 1 


21 


x^ + x^ - 6 x^ - 6 x^ + 8 x^ + 8 X + 1 


22 


x^*^ - 11 X** + 44 x*^ - 77 x^ + 55 x^ - 11 


23 


x" - x^" - 10 x^ + 9 x« + 36 x' - 28 x^ - 56 x^ + 35 x'^+ 
35x^ - 15x^ - 6x + 1, 


24 


X* - 8x'^ + 20x'' - 16x^ + 1 




^lU 1 fi ^8 _|_ qcr ™6 ™5 ™.4 i cr 'S i ™2 ft ™ i 


26 


x^^ - 13x1" + 65x« - 156x^ + 182x'' - 91 x^ + 13 


27 


x'-* - 9x'^ + 27 x^ - 30x^ + 9x - 1 


28 


xi2 - 12x1° ^ 53^8 _ ;lo4x6 + 86x4 - 24^x'^ + 1 


29 


^14 _ ^13 _ I'^^vi 12 x^^ + 66x1" _ 55^9 _ 165 x«+ 
120x^ + 210x^ - 126x5 - 126x4 + 56x3 + 28x2 - 7x - 1 


30 


X* - 7x'^ + 14x4 -8x^ + 1 







Table 3: A187360 fn, m) coefficient array of 
minimal polynomials of 2 cos ^— j , rising powers 



n/m 





1 


2 


3 


4 


5 


6 




1 


2 


1 














2 





1 














3 


-1 


1 














4 


-2 





1 












5 


-1 


-1 


1 












6 


-3 





1 












7 


1 


-2 


-1 


1 










8 


2 





-4 





1 








9 


-1 


-3 





1 










10 


5 





-5 





1 








11 


-1 


3 


3 


-4 


-1 


1 






12 


1 





-4 





1 








13 


-1 


-3 


6 


4 


-5 


-1 


1 




14 


-7 





14 





-7 





1 




15 


1 


-4 


-4 


1 


1 
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Table 4: Zeros of C(n, x) in power basis (rising powers of p(n)) for n = 1,2,. ..,30. 



n 


coemcients or C-zeros in power basis < p , ' > 


1 


[[-2]] 


2 


[[0]] 


3 


[[1]] 


4 


[[0, i],[o,-i]] 


5 


[[0, 


6 


[[0, i],[o,-i]] 


7 


[[0, 1], [-1,-1,1], [2, 0,-1]] 


8 


[[0, 1], [0,-3, 0,1], [0,3, 0,-1], [0,-1]] 


9 


[[0, 1], [-2, -1,1], [2, 0,-1]] 


10 


[[0, 1], [0,-3, 0,1], [0,3, 0,-1], [0,-1]] 


11 


[[0, 1], [0, -3, 0, 1], [1,2, -3, -1, 1], [-2, 0, 4, 0,-1], [2, 0, -1]] 


12 


[[0, 1], [0,4, 0,-1], [0,-4,0, 1],[0,-1]] 


13 


[[0, 1], [0, -3, 0, 1], [0, 5, 0, -5, 0, 1], [1, -3, -3, 4, 1, -1], [-2, 0, 4, 0, -1], [2, 0, -1]] 


14 


[[0, 1], [0, -3, 0, 1], [0, 5, 0, -5, 0, 1], [0, -5, 0, 5, 0, -1], [0, 3, 0, -1], [0, -1]] 


15 


[[0, 1], [-2, 3, 1,-1], [-1,-4, 0,1], [2, 0,-1]] 


16 


[[0, 1], [0, -3, 0, 1], [0, 5, 0, -5, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [0, 7, 0, -14, 0, 7, 0, -1], [0, -5, 0, 5, 0, -1], 
[0,3, 0,-1], [0,-1]] 


17 


[[0, 1], [0, -3, 0, 1], [0, 5, 0, -5, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [-1,4, 6, -10, -5, 6, 1, -1], 
[2, 0, -9, 0, 6, 0, -1], [-2, 0, 4, 0, -1], [2, 0, -1]] 


18 


[[0, 1], [0, 5, 0, -5, 0, 1], [0, -4, 0, 5, 0,-1], [0, 4, 0, -5, 0, 1], [0, -5, 0, 5, 0,-1], [0, -1]] 


19 


[[0, 1], [0, -3, 0, 1], [0, 5, 0, -5, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [1, 4, -10, -10, 15, 6, -7, -1, 1], 
[-2, 0, 16, 0, -20, 0, 8, 0, -1], [2, 0, -9, 0, 6, 0, -1], [-2, 0, 4, 0, -1], [2, 0, -1]] 


20 


[[0, 1], [0, -3, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [0, 8, 0, -18, 0, 8, 0, -1], [0, -8, 0, 18, 0, -8, 0, 1], 
[0, 7, 0, -14, 0, 7, 0, -1], [0, 3, 0, -1], [0, -1]] 


21 


[[0, 1], [0, 5, 0, -5, 0, 1], [2, 3, -4, -1, 1], [-3, -9, 1, 6, 0, -1], [-2, 0, 4, 0,-1], [2, 0, -1]] 


22 


[[0, 1], [0, -3, 0, 1], [0, 5, 0, -5, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [0, 9, 0, -30, 0, 27, 0, -9, 0, 1], 

[0, -9, 0, 30, 0, -27, 0, 9, 0,-1], [0, 7, 0, -14, 0, 7, 0,-1], [0, -5, 0, 5, 0,-1], [0, 3, 0,-1], [0, -1]] 


23 


[[0, 1], [0, -3, 0, 1], [0, 5, 0, -5, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [0, 9, 0, -30, 0, 27, 0, -9, 0, 1], 
[-1, -5, 15, 20, -35, -21, 28, 8, -9, -1, 1], [2, 0, -25, 0, 50, 0, -35, 0, 10, 0,-1], 
[-2, 0, 16, 0, -20, 0, 8, 0, -1], [2, 0, -9, 0, 6, 0, -1], [-2, 0, 4, 0, -1], [2, 0, -1]] 


24 


[[0, 1], [0, 5, 0, -5, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [0, -8, 0, 6, 0, -1], [0, 8, 0, -6, 0, 1], 
[0, 7, 0, -14, 0, 7, 0, -1], [0, -5, 0, 5, 0, -1], [0, -1]] 


25 


[[0, 1], [0, -3, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [0, 9, 0, -30, 0, 27, 0, -9, 0, 1], 

[n —10 ^ '^n — ^ —97 1 Q — ll [n in — 1^ — 1^ 90 7 —R — l ll [—9 Ifi —90 ^0 — ll 
[2, 0, -9, 0, 6, 0, -1], [-2, 0, 4, 0, -1], [2, 0, -1]] 


26 


[[0, 1], [0, -3, 0, 1], [0, 5, 0, -5, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [0, 9, 0, -30, 0, 27, 0, -9, 0, 1], 

[0, -11, 0, 55, 0, -77, 0, 44, 0, -11,0, 1], [0, 11, 0, -55, 0, 77, 0, -44, 0, 11, 0, -1], 

[0, -9, 0, 30, 0, -27, 0, 9, 0, -1], [0, 7, 0, -14, 0, 7, 0, -1], [0, -5, 0, 5, 0, -1], [0, 3, 0, -1], [0, -1]] 


27 


[[0, 1], [0, 5, 0, -5, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [-2, 7, 1, -14, 0, 7, 0, -1], [2, -5, -4, 5, 1, -1], 
[2, -1, -16, 0, 20, 0, -8, 0, 1], [-2, 0, 16, 0, -20, 0, 8, 0, -1], [-2, 0, 4, 0, -1], [2, 0, -1]] 




continued on next page 



Note: the higher power coefficients not shown are all zero. 

Example: n = 27: the first zero of C(27,x) is p(27), the second zero is 5p(27) — 5 p(27)^ — 
lp(27)^, etc. 
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Table 4 continued: Zeros of C(n, x) in power basis (rising powers of p(n)) for 

n = 28, ...,30. 



n 


power basis coefficients for < p^, p''^") ^ > 


28 


[[0, 1], [0, -3, 0, 1], [0, 5, 0, -5, 0, 1], [0, 9, 0, -30, 0, 27, 0, -9, 0, 1], [0, -11, 0, 55, 0, -77, 0, 44, 0, -11,0, 1], 
0, 12, 0, -67, 0, 96, 0, -52, 0, 12, 0, -1], [0, -12, 0, 67, 0, -96, 0, 52, 0, -12, 0, 11, 

[0, 11,0, -55, 0, 77, 0, -44, 0, 11, 0, -1], [0, -9, 0, 30, 0, -27, 0, 9, 0, -1], [0, -5, 0, 5, 0, -1], [0, 3, 0, -1], 

[0,-1]] 


29 


[[0, 1], [0, -3, 0, 1], [0, 5, 0, -5, 0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [0, 9, 0, -30, 0, 27, 0, -9, 0, 1], 

[0, -11, 0, 55, 0, -77, 0, 44, 0, -11,0, 1], [0, 13, 0, -91, 0, 182, 0, -156, 0, 65, 0, -13, 0, 1], 

[1, -7, -21, 56, 70, -126, -84, 120, 45, -55, -11, 12, 1, -1], [-2, 0, 36, 0, -105, 0, 112, 0, -54, 0, 12, 0,-1], 

[2, 0, -25, 0, 50, 0, -35, 0, 10, 0, -1], [-2, 0, 16, 0, -20, 0, 8, 0, -1], [2, 0, -9, 0, 6, 0, -1], [-2, 0, 4, 0, -1], 

[2,0,-1]] 


30 


[[0, 1], [0, -7, 0, 14, 0, -7, 0, 1], [0, -7, 0, 22, 0, -13, 0, 2], [0, 4, 0, -13, 0, 7, 0,-1], [0, -4, 0, 13, 0, -7, 0, 1], 
[0, 7, 0, -22, 0, 13, 0, -2], [0, 7, 0, -14, 0, 7, 0, -1], [0, -1]] 
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Table 5: Extended set Ai(n), first differences AAl(n), and floor-pattern ^(n) for composed odd squarefree 

modulus n. 



m 


n(m) 


2S{n) 


M{n) 


AM{n) 


^(n) 


1 


15 


8 


[0,1,7,11,13,17] 


[1,6,4,2,4] 


[2,1,0,1,0,1,2,0] 


2 


21 


12 


[0,1,5,11,13,17,19,23] 


[1,4,6,2,4,2,4] 


[1,2,0,1,0,1,0,1,0,2,1,0] 


3 


33 


20 


[0, 1, 5, 7, 13, 17, 19, 23, 25, 29, 31, 35] 


[1,4,2,6,4,2,4,2,4,2,4] 


[1,0,2,1,0,1,0,1,0,1, 
0,1,0,1,0,1,2,0,1,0] 


4 


35 

5 • 7 


24 


[0, 1, 3, 9, 11, 13, 17, 19, 23, 27, 29, 31, 33, 37] 


[1,2,6,2,2,4,2,4,4,2,2,2,4] 


[0,2,0,0,1,0,1,1,0,0,0,1, 
0,0,0,1,1,0,1,0,0,2,0,0] 


5 


39 

3 • 13 


24 


[0, 1, 5, 7, 11, 17, 19, 23, 25, 29, 31, 35, 37, 41] 


[1,4,2,4,6,2,4,2,4,2,4,2,4] 


[1,0,1,2,0,1,0,1,0,1,0,1, 
0,1,0,1,0,1,0,2,1,0,1,0] 


6 


51 

3 • 17 


32 


[0,1,5,7,11,13,19,23,25,29, 
31,35,37,41,43,47,49,53] 


[1,4,2,4,2,6,4,2,4,2,4,2,4,2,4,2,4] 


[1,0,1,0,2,1,0,1,0,1,0,1,0,1,0,1, 
0,1,0,1,0,1,0,1,0,1,2,0,1,0,1,0] 


7 


55 
5 • 11 


40 


[0, 1, 3, 7, 9, 13, 17, 19, 21, 23, 27, 29, 
29, 31, 37, 39, 41, 43, 47, 49, 51, 53, 57] 


[1,2,4,2,4,4,2,2,2,4,2, 
2,6,2,2,2,4,2,2,2,4] 


[0,1,0,1,1,0,0,0,1,0,0,2,0,0, 
0,1,0,0,0,1,0,0,0,1,0,0,0,2,0, 
0,1,0,0,0,1,1,0,1,0,0] 


8 


57 

3 • 19 


36 


[0,1,5,7,11,13,17,23,25,29,31, 
31, 35, 37, 41, 43, 47, 49, 53, 55, 59] 


[1,4,2,4,2,4,6,2,4, 
2,4,2,4,2,4,2,4,2,4] 


[1,0,1,0,1,2,0,1,0,1,0,1,0,1,0, 
1,0,1,0,1,0,1,0,1,0,1,0,1,0, 
2,1,0,1,0,1,0] 


9 


65 

5 • 13 


48 


[0, 1, 3, 7, 9, 11, 17, 19, 21, 23, 27, 29, 31, 33, 
33, 37, 41, 43, 47, 49, 51, 53, 57, 59, 61, 63, 67] 


[1,2,4,2,2,6,2,2,2,4,2,2, 
2,4,4,2,4,2,2,2,4,2,2,2,4] 


[0,1,0,0,2,0,0,0,1,0,0,0,1,1,0,1, 
0,0,0,1,0,0,0,1,0,0,0,1,0,0,0,1, 
0,1,1,0,0,0,1,0,0,0,2,0,0,1,0,0] 


10 


69 

3 • 23 


44 


[0, 1, 5, 7, 11, 13, 17, 19, 25, 29, 31, 35, 37, 
41, 43, 47, 49, 53, 55, 59, 61, 65, 67, 71] 


[1,4,2,4,2,4,2,6,4,2,4,2, 
4,2,4,2,4,2,4,2,4,2,4] 


[1,0,1,0,1,0,2,1,0,1,0,1,0,1,0, 
1,0,1,0,1,0,1,0,1,0,1,0,1,0,1, 
0,1,0,1,0,1,2,0,1,0,1,0,1,0] 


11 


77 

7 • 11 


60 


[0, 1, 3, 5, 9, 13, 15, 17, 19, 23, 25, 27, 
29, 31, 37, 39, 41, 43, 45, 47, 51, 53, 
57, 59, 61, 65, 67, 69, 71, 73, 75, 79] 


[1,2,2,4,4,2,2,2,4,2,2,2,2,6,2,2, 
2,2,2,4,2,4,2,2,4,2,2,2,2,2,4] 


[0,0,1,1,0,0,0,1,0,0,0,0,2,0,0, 
0,0,0,1,0,1,0,0,1,0,0,0,0,0,1, 
0,0,0,0,0,1,0,0,1,0,1,0,0,0,0, 
0,2,0,0,0,0,1,0,0,0,1,1,0,0,0] 


continued 





Table 5 ctnd.: Extended set Al(n), first differences A A4(n), and fioor-pattern ^(n) for composed odd squarefree 

modulus n. 



in 


n(m) 


2S{n) 


Min) 


AM{n) 




12 
13 


85 

5 • 17 

87 

3 • 29 


64 
56 


[0, 1, 3, 7, 9, 11, 13, 19, 21, 23, 27, 29, 31, 
33, 37, 39, 41, 43, 47, 49, 53, 57, 59, 61, 
63, 67, 69, 71, 73, 77, 79, 81, 83, 87] 

[0,1,5,7,11,13,17,19,23,25,31 
35, 37, 41, 43, 47, 49, 53, 55, 59, 61, 
65,67,71,73,77, 79,83,85,89] 


[1,2,4,2,2,2,6,2,2,4,2, 
2,2,4,2,2,2,4,2,4,4,2, 
2,2,4,2,2,2,4,2,2,2,4] 

[1,4,2,4,2,4,2,4,2,6,4,2,4,2,4, 
2,4,2,4,2,4,2,4,2,4,2,4,2,4] 


[0,1,0,0,0,2,0,0,1,0,0,0,1,0,0,0, 
1,0,1,1,0,0,0,1,0,0,0,1,0,0,0,1, 
0,0,0,1,0,0,0,1,0,0,0,1,1,0,1,0, 
0,0,1,0,0,0,1,0,0,2,0,0,0,1,0,0] 
[1,0,1,0,1,0,1,0,2,1,0,1,0,1, 
0,1,0,1,0,1,0,1,0,1,0,1,0,1 
0,1,0,1,0,1,0,1,0,1,0,1,0,1, 
0,1,0,1,2,0,1,0,1,0,1,0,1,0] 


14 


91 

7 • 13 


72 


[0,1,3,5,9,11,15,17,19,23,25, 
27, 29,31,33,37,41,43,45,47, 
51,53,55,57, 59,61,67,69,71, 
73,75,79,81,83,85,87,89,93] 


[1,2,2,4,2,4,2,2,4,2,2,2,2, 
2,4,4,2,2,2,4,2,2,2,2,2,6, 
2,2,2,2,4,2,2,2,2,2,4] 


[0,0,1,0,1,0,0,1,0,0,0,0,0,1,1, 
0,0,0,1,0,0,0,0,0,2,0,0,0,0,1, 
0,0,0,0,0,1,0,0,0,0,0,1,0,0,0, 
0,0,0,1,0,0,1,0,1,0,0,0] 


15 


93 

3 • 31 


60 


[0, 1, 5, 7, 11, 13, 17, 19, 23, 25, 29, 35, 
37, 41, 43, 47, 49, 53, 55, 59, 61, 65, 67, 
71,73,77, 79,83,85,89,91,95] 


[1,4,2,4,2,4,2,4,2,4,6, 
2,4,2,4,2,4,2,4,2,4,2, 
4,2,4,2,4,2,4,2,4] 


[1,0,1,0,1,0,1,0,1,2,0,1,0,1,0, 
1,0,1,0,1,0,1,0,1,0,1,0,1,0,1, 
0,1,0,1,0,1,0,1,0,1,0,1,0,1,0, 
1,0,1,0,2,1,0,1,0,1,0,1,0,1,0] 


16 


95 

5 • 19 


72 


[0,1,3,7,9,11,13,17,21,23,27, 
29,31,33,37, 39,41,43,47, 49, 

51,53,59,61,63,67,69,71,73, 

77 70 ai a"? a7 so qi q"? q71 
/ ( , /y, oi, oo, o ( , oy, yi, yo, y / j 


[1,2,4,2,2,2,4,4,2,4,2,2,2, 
4,2,2,2,4,2,2,2,6,2,2,4,2, 
2,2,4,2,2,2,4,2,2,2,4] 


[0,1,0,0,0,1,1,0,1,0,0,0,1,0,0, 
0,1,0,0,0,2,0,0,1,0,0,0,1,0,0, 
0,1,0,0,0,1,0,0,0,1,0,0,0,1,0, 

U,U, i,U,U,Z,U,U,U, i,U,U,U, i,u, 

0,0,1,0,1,1,0,0,0,1,0,0] 


17 


105 

3 • 5 • 7 


48 


[0,1,11,13,17,19,23,29,31,37, 
41,43,47, 53,59,61,67,71,73, 
79,83,89,97, 101,103,107] 


[1,10,2,4,2,4,6,2,6, 
4,2,4,6,6,2,6,4,2, 
6,4,6,8,4,2,4] 


[4,0,1,0,1,2,0,2,1,0,1,2, 
2,0,2,1,0,2,1,2,3,1,0,1, 
0,1,3,2,1,2,0,1,2,0,2,2, 
1,0,1,2,0,2,1,0,1,0,4,0] 















Table 6: Cycle structure of ^a/(Q(p(n))/Q) for n = 1,2, ...,40. 



n 


cycles 


1 


[[0]] = [[1]] 

LLJJ LLJJ 


2 


[[1]] 

L L J J 


3 


[[1]] 


4 


[[3, 1]] 


5 


[[3, 1]] 


6 


[[5,1]] 


7 


[[3,5,1]] 


8 


[[3,7,5,1]] 


9 


[[5,7, 1]] 


10 


[[3,9,7, 1]] 


11 


[[3,9,5,7, 1]] 




12 




[[5,1], [7, 1],[11,1]] 


13 


[[7,3,5,9,11,1]] 


14 


[[5,3,13,9,11,1]] 


15 


[[7,11,13, 1]] 


16 


[[3,9,5,15,13,7,11,1]] 


17 


[[3,9,7,13,5,15,11,1]] 


18 


[[5,11,17,13,7, 1]] 


19 


[[3,9,11,5,15,7,17, 13,1]] 




20 




[[3,9,13, 1],[7,9,17, 1], [11,1], [19,1]] 


21 


[[11,5,13,17, 19,1]] 


22 


[[3,9,17,7,21,19,13,5,15,1]] 


23 


[[3,9,19,11,13,7,21,17,5,15,1]] 




24 




[[5,23,19, 1],[7, 1], [11, 23, 13,1], [17,1]] 


25 


[[3,9,23,19,7,21,13,11,17, 1]] 


26 


[[7, 3, 21, 9, 11, 25, 19, 23, 5, 17, 15, 1]] 


27 


[[5,25,17, 23,7,19,13,11,1]] 




28 




[[3, 9, 27, 25, 19, 1], [5, 25, 13, 9, 11, 1], [17, 9, 15, 25, 23, 1]] 


29 


[[3, 9, 27, 23, 11, 25, 17, 7, 21, 5, 15, 13, 19, 1]] 




30 




[[7, 11, 17, 1], [13, 11, 23, 1], [19, 1], [29, 1]] 


31 


[[3, 9, 27, 19, 5, 15, 17, 11, 29, 25, 13, 23, 7, 21, 1]] 


32 


[[3, 9, 27, 17, 13, 25, 11, 31, 29, 23, 5, 15, 19, 7, 21, 1]] 


33 


[[5,25,7,31,23,17, 19,29,13, 1]] 


34 


[[3, 9, 27, 13, 29, 19, 11, 33, 31, 25, 7, 21, 5, 15, 23, 1]] 


35 


[[3, 9, 27, 11, 33, 29, 17, 19, 13, 31, 23, 1]] 




36 




[[5, 25, 19, 23, 29, 1], [7, 23, 17, 25, 31, 1], [11, 23, 35, 25, 13, 1]] 


37 


[[5, 25, 23, 33, 17, 11, 19, 21, 31 , 7, 35, 27, 13, 9, 29, 3, 15, 1]] 


38 


[[13, 17, 7, 15, 33, 27, 29, 3, 37, 25, 21, 31, 23, 5, 11, 9, 35, 1]] 


39 


[[7, 29, 31, 17, 37, 25, 19, 23, 5, 35, 11, 1]] 




40 




[[3,9,27, 1], [7,31,23, 1], [11,39,29, 1], [13,9,37, 1], [17,31,33, 1], [19,39,21,1]] 







Boxed and colored n-numbers indicate non-cyclic Galois groups. See Table 7. 
Boldface numbers are nontrivial square roots Modd n, denoted by s in the text. 
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Table 7: Non-cyclic Galois groups ^a/(Q(C(n))/Q) , n < 100 



n 


ip{n) 


cycle structure 


no. of cycles 


Galoisgroup 




8 


4 


23 


3 


Z2 X Z2 


12 


4 


23 


3 


Z2 X Z2 




15 




8 


h 22 


4 


Z4 X Z2 


16 


8 


42 22 


4 


Z4 X Z2 


20 


8 


42 22 


4 


Z4 X Z2 




21 




12 


63 


3 


Z3 X Z| 




24 




8 


27 


7 


Z| 


28 


12 


63 


3 


Z3 X 


30 


8 


42 22 


4 


Z4 X Z2 




32 




16 


82 4i 22 


5 


Zs X Z2 




33 




20 


103 


3 


Z5 X Z| 




35 




24 


122 62 


4 


Z4 X Z3 X Z2 


36 


12 


63 


3 


Z3 X Z'i 


39 


24 


122 62 


4 


Z4 X Z3 X Z2 




40 




16 


44 26 


10 


Z4 X zl 


42 


12 


63 


3 




44 


20 


103 


3 


Z5 X 


45 


24 


122 62 


4 


Z4 X Z3 X Z2 


48 


16 


44 26 


10 


Z4 X Z| 




51 




32 


162 8i4i 22 


6 


Z16 X Z2 


52 


24 


122 62 


4 


Z4 X Z3 X Z2 




55 




40 


202 102 


4 


Z5 X Z4 X Z2 




56 




24 


67 


7 


Z3 X Z| 




57 




36 


183 


3 


Zg X Z| 


60 


16 


44 26 


10 


Z4 X Z'l 




63 




36 


612 


12 


X 


64 


32 


162 8i4i 22 


6 


Z16 X Z2 




65 




48 


126 


6 


Z| X Z3 


66 


20 


IO3 


3 


Z5 X Z'l 


68 


32 


162 8i4i 22 


6 


Z16 X Z2 




69 




44 


223 


3 


Zii X Z| 


70 


24 


122 62 


4 


Z4 X Z3 X Z2 


72 


24 


67 


7 


Z3 X Z| 


75 


40 


2O2 IO2 


4 


Z5 X Z4 X Z2 


76 


36 


I83 


3 


Zg X Z| 




77 




60 


3O3 


3 


Z5 X Z3 X Z| 


78 


24 


122 62 


4 


Z4 X Z3 X Z2 




80 




32 


4i2 24 


16 


X Z2 



Continued on the next page. 
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Table 7 continued: Non-cyclic Galois groups ^a/(Q(<^(n))/Q) , n < 100 



n 


d{n) 


cycle structure 


no. of cycles 


Galoisgroup 


84 


24 


67 


7 


Zs X Zl 




85 




64 


I64 82 44 


10 


ZlQ X Z4 




87 




56 


282 142 


4 


Zj X Z4 X Z2 




88 




40 


107 


7 


Z5 X Z| 


90 


24 


122 62 


4 


Z4 X Z3 X Z2 




91 




72 


128 68 


16 


Zi X zl X Z2I 




92 




44 


223 


3 


r-r r-rO 

Zii X Zl 


93 


60 


303 


3 


Zf, X Zo,x Zl 




95 




72 


362 182 


4 


Zq X Z4 X Z2 




96 




32 


8443 26 


13 


Zs X Zl 


99 


60 


303 


3 


Z5 X Z3 X zl 


100 


40 


202 102 


4 


Z^ X Z4 X Z2 













The cyclic group of order m is denoted by Zm. For all other values n < 100 the Galois group is the 
cyclic group Z^(n) . 

Only independent cycles are counted, i.e., cycles which appear as sub-cycles of the given ones have been 
omitted. 

The notation, e.g., I62 81 4i 22, means that there are 2 cycles of order (length) 16, one cycle of order 8, 
one cycle of order 4 and two cycles of order 2. 

Direct products of identical cyclic groups are sometimes written in exponent form, e.g., Z2 stands for 
Z2 X Z2. 

Boxed and colored n-numbers indicate where some non-cyclic Galois group appears for the first time. 
Some of the cycle graphs are shown in Fig. 4. 



54 



Table 8: Non-cyclic Galois groups ^a/(Q(p(n))/Q) , n < 100 



n 


Sin) 


cycle structure 


no. of cycles 


Galoisgroup 




12 




4 


23 


3 


Z2 X Z2 




20 




8 


h 22 


4 


Z4 X Z2 


24 


8 


h 22 


4 


Z,i X Z2 




28 




12 


63 


3 


Z3 X Zi 


30 


8 


42 22 


4 


Zi X Z2 


36 


12 


63 


3 


Z3 X Z'i 




40 




16 


46 


6 


Zi X Z4 


42 


12 


63 


3 


Z3 X Zl 




44 




20 


103 


3 


Z5 X Z| 




48 




16 


82 4i 22 


5 


Zg X Z2 




52 




24 


122 62 


4 


Z4 X Z3 X Z2 


56 


24 


122 62 


4 


Z4 X Z3 X Z2 




60 




16 


44 26 


10 


Z4 X Z| 




63 




18 


64 


4 


Zl X Z2 


65 


24 


122 62 


4 


Z4 X Z3 X Z2 


66 


20 


103 


3 


Z5 X 




68 




32 


I62 81 4i 22 


6 


Z16 X Z2 


70 


24 


122 62 


4 


Z4 X Z3 X Z2 


72 


24 


122 62 


4 


Z4 X Z3 X Z2 




76 




36 


183 


3 


Zg X Zl 


78 


24 


122 62 


4 


Z4 X Z3 X Z2 




80 




32 


84 44 


8 


Zs X Z4 




84 




24 


67 


7 


Z3 X zl 


85 


32 


I62 81 4i 22 


6 


Z16 X Z2 




88 




40 


202 102 


4 


Z5 X Z4, X Z2 


90 


24 


122 62 


4 


Z4 X Z3 X Z2 




91 




36 


124 


4 


Zi X zl 


92 


44 


223 


3 


Zii X Zl 


96 


32 


I62 81 4i 22 


6 


Z16 X Z2 


100 


40 


2O2 IO2 


4 


Z5 X Z4 X Z2 













The cyclic group of order m is denoted by Zm- For all other n < 100 cases the Galois group is the 
cyclic group Z^^^). The n values are given in A206552 . 

Only independent cycles are counted, i.e., cycles which appear as sub-cycles of the given ones have been 
omitted. 

The notation, e.g., I62 81 4i 22, means that there are 2 cycles of order (length) 16, one cycle of order 8, 
one cycle of order 4 and two cycles of order 2. 

Direct products of identical cyclic groups are sometimes written in exponent form, e.g., Z2 stands for 
Z2 X Z2. 

Boxed and colored n-numbers indicate where some non-cyclic Galois group appears for the first time. 
For the cycle graphs see Fig. 4. 
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8=4: Z_2 X Z_2 



8=8: Z_4xZ_2 



8=12: Z_3 X Z_2 x Z_2 



8=16: Z_4xZ_4 







8=20: Z_5 X Z_2 x Z_2 



8=16:Z_8xZ_2 



8=24: Z_4 x Z_3 x Z_2 



8=16: Z_4 x Z_2 x Z_2 






8=18: Z 3xZ 3xZ 2 



8=32: Z 16xZ 2 



8=36: Z 9xZ 2xZ 2 



Figure 4: Cycle graphs for non-cyclic Galois groups Q al{'Q{p{n)) /Q) appearing for n = 1..100. 5 is the 
degree of C(n,x), the maximal polynomial of p{n), hence the order of the Galois group. See Table 8. 
Continued on next page. 
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6=32: Z 8 xZ 4 



5=24: Z 3xZ 2xZ 2xZ 2 5=40: Z 5 x Z 4 x Z 2 




• • • 



5=36: Z4xZ3xZ3 



Figure 4 continued: Cycle graphs for non-cyclic Galois groups Gal{Q{p{n 
See Table 8. 



appearing for n = 1..100. 
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